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Abstract 



Gauge/Gravity dualities open the non-perturbative realms of strongly-coupled 
gauge theories to analytic treatment. Anti-de Sitter Space/Conformal Field The- 
ory, one way of connecting gravity dual models to gauge theories, is a correspon- 
dence between a ten-dimensional Type IIB superstring theory in AdS^ x and a 
A/" = 4 super Yang Mills theory. To describe systems that are experimentally ac- 
cessible, however, the formal correspondence is modified into a phenomenological 
duality between a five-dimensional gravity model and a strongly coupled QCD-like 
gauge theory. This duality is referred to as AdS/QCD. 

This work explores aspects of the soft-wall AdS/QCD model. The phrase 'soft 
wall' refers to the means of breaking the conformal symmetry and introducing 
a mass scale to the gauge side of the duality. We add higher-order terms to 
the soft-wall Lagrangian and calculate the effect on physical observables. Meson 
mass spectra gain a more complex structure, exhibiting a better match with the 
experimental values than previous models. The Gell-Mann-Oakes-Renner rela- 
tion naturally emerges from the model. We calculate the form factor and the 
coupling QpT^T^ as a non-trivial test on the limits of our soft-wall model. 

Introducing a black brane into the gravity dual metric allows us to derive 
thermodynamic quantities in the gauge theory. As expected at high temperatures, 
the entropy scales as the cube of the temperature, and the speed of sound reaches 
its conformal limit of one-third. Thermal condensates contribute leading-order 
terms, modifying the temperature, entropy, and free energy behavior. We find 
that the system undergoes a phase transition from hadronic matter to a strongly 
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coupled quark-gluon plasma at a critical temperature. 
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Chapter 1 
Introduction 



"You have to have an idea of what you are going to do, but it should be a vague 

idea." 
-Pablo Picasso 

All the matter and energy in the 100 billion stars of the Milky Way galaxy 
and the 100 billion galaxies in the universe came from one infinitely dense point. 
Some 13.1 to 13.7 billion years ago, the point began expanding after the event 
now known as the Big Bang . According to modern thinking, one unified force 
governed the universe for the first 10"'^^ seconds of its existence [HI [1]. Over the 
next picosecond (10^^^ s), the unified force fractured into the four known forces 
that rule the cosmos today: strong, electromagnetic, weak, and gravitational. 

Today, we have two prevailing theories describing all of nature. The Standard 
Model encapsulates the strong, electromagnetic, and the weak forces, describing 
them in terms of so-called gauge theories [12]. General Relativity describes the 
spacetime curvature responsible for gravitational interactions, but is inherently 
incompatible with the Standard Model. Many physicists, past and present, seek 
theoretical constructs that unify the Standard Model with general relativity into 
a grand unified theory that depicts the beginning of the universe. One such grand 
unified scheme that has gained notoriety for its mathematical complexity and ex- 
clusion from experimental verification is string theory. However, ideas from string 
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theory allow for a correspondence between extra-dimensional gravity models and 
strongly coupled gauge theories. Such a correspondence allows calculations done 
in gravity models to be interpreted as gauge theory quantities. In our research, we 
concentrate on how gauge/gravity dualities describe the strongly coupled regions 
of the strong force. 

A gauge theory called quantum chromodynamics (QCD) describes the strong 
force. QCD characterizes the interactions among quarks, the constituents of nu- 
cleons, and gluons, the color-charge force carriers. Using QCD, one can develop 
effective theories describing hadronic matter. QCD distinguishes itself from other 
gauge theories with two characteristics: confinement and asymptotic freedom. 
Confinement ensures that only color-charge singlets exist, meaning only color- 
less objects are physically allowed. Asymptotic freedom refers to the high-energy 
behavior of quarks. At high energies, the quarks within the hadrons are weakly- 
interacting particles, allowing perturbative calculations. 

Traditionally, calculations within gauge theories have been done in series ex- 
pansions. For instance, if we wanted to calculate the cross section of an interaction 
CT/, then we could express it as a perturbative expansion in powers of the coupling 
constant g, where q is the energy-momentum vector of the collision, 

aj = F,{q')g + F2{q')g^ + F,{q')g' + .... (1.1) 

As long as g <^ 1, higher-order terms in g may be ignored, and calculations 
consider only the first few terms. A problem arises when g > I; subsequent 



terms add ever-increasing values to the expansion (1.1 ). Thus for strongly coupled 
systems {g > 1), perturbative techniques break down. 

In QCD, as in all gauge theories in the Standard Model, the coupling con- 
stant changes with energy, g — )■ g{q^); this phenomenon is termed the running of 
the coupling. As the energies involved increase, the coupling constant of QCD 
decreases. At energies less than 0{1 GeV), perturbative calculations fail, and 
new means of calculating physical observables must be found [13] . Gauge/gravity 
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correspondence provides a remarkable mathematical tool to carry out such non- 
perturbative calculations in strongly coupled theories. 
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Figure 1.1: The phase diagram of QCD [Ij. Hadronic matter experiences a cross- 
over transition to quark-gluon plasma as it reaches a critical temperature. Exper- 
iments at the RHIC and LHC probe temperatures around the cross-over. 

According to recent experiments, nonperturbative calculations are still needed 
as the temperature is increased in physical systems. In QCD, the protons and neu- 
trons melt and become quark-gluon plasma (QGP) with increasing temperatures. 
Just as water turns to steam, as temperatures increase beyond a critical point, 
the confined hadronic matter undergoes a cross-over phase change to deconfined 
state of the QGP as seen in the phase diagram in Figure [TTT| In this state, quarks 
and gluons are no longer bound within hadrons but associate freely [HllTn]. Such 
a substance has been created at the Relativistic Heavy-Ion Collider (RHIC) at 
Brookhaven National Laboratory and most recently at the Large Hadron Collider 
(LHC) at CERN. Experimental results suggest that the QGP is strongly coupled 
and mimics the behavior of fluids [lU HTJ dSl [12]. As a result, ideas from fluid 
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dynamics are useful in gauge/gravity dualities at finite temperature. 

In this thesis, we will discuss a particular gauge/gravity duality called Anti-de 
Sitter /Conformal Field Theory (AdS/CFT) and how it gives rise to an effective 
correspondence called Anti-de Sitter/Quantum Chromodynamics (AdS/QCD). 
This work proceeds as follows 

• In Chapter [2] and [3} we review the literature concerning the correspondence 
dictionary of AdS/CFT and explain how it yields the effective correspon- 
dence AdS/QCD. We go into detail on how to modify the gravity dual model 
in AdS space to describe a QCD-like theory with a large number of colors, 
Nc- We explore how physicists have applied the gauge/gravity duality to 
model condensates, chiral symmetry breaking (CSB), and meson mass spec- 
tra of QCD. We also detail how the duality characterizes the QGP and 
how the thermodynamics of this QGP-like fluid are encoded in a black-hole 
metric of the gravity model. 

• In Chapter |4| we modify existing AdS/QCD models at zero temperature. 
We improve the description of CSB done in previous works. These improve- 
ments are expanded to the mass spectra of the scalar, pseudoscalars, vectors, 
and the axial- vectors. Concluding this chapter, we explore AdS/QCD cal- 
culations of QpTTTT, the form factor F^^, and the Gell-Mann-Oakes-Renner 
relation. 

• In Chapter [sj we expand our AdS/QCD model by introducing a black brane 
into the AdS space. We then explore the entropy, speed of sound, and 
free energy of the QGP-like fluid created. Investigating the actions of the 
systems, we show how the transition from the AdS thermal metric to a 
black-hole metric corresponds to the transition from the confined phase to 
the deconfined phase of QCD. 

• Finally, in Chapter [6} we discuss the main conclusions of the analysis, open 
questions yet to be answered, and future prospects for this area of research. 



The majority of the work presented here has been previously pubhshed in 
and [21j. The material in Chapter [s] and Appendix C.3 will be published at a 
later date. Additional details and expanded discussions are included. 



Chapter 2 

AdS/CFT: Overview of the 
Gauge/Gravity Duality 

"This is your last chance. After this, there is no turning back. You take the blue 
pill - the story ends, you wake up in your bed and believe whatever you want to 
believe. You take the red pill - you stay in Wonderland and I show you how deep 

the rabbit-hole goes." 
Morpheus, The Matrix{1999) 

In this chapter, we explain how the AdS/CFT correspondence connects to 
QCD. The focus of this thesis is modeling particular dynamics and thermody- 
namics of a Q CD-like gauge theory using a soft-wall AdS/QCD. Explaining the 
connection between AdS / CFT to AdS / QCD and then the phenomenology leading 
up to the original research in Chapter |4] and [5] covers a great deal of material. For 
additional background, several reviews have been written [221 ESI El 1211 12SI [2S] • 
We begin by exploring the basics of the AdS/CFT correspondence. From there, 
we introduce several phenomenological modifications meant to deform the CFT 
into a QCD-like theory with a large number of colors. 
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2.1 The Basics of AdS/CFT 



The basis of all dualities is the ability to describe a physical system in two dif- 
ferent but equivalent ways. The principles of AdS/CFT correspondence were 
established late last century [271 ESI ESI ED]- The seminal works drew a direct 
correspondence between Type JIB string theory in AdS^ x with a A/" = 4 su- 
persymmetric Yang-Mills (SYM) theory by investigating the behavior of so-called 
branes in the low-energy and strong coupling limits. Because AdS/CFT uses 
a lower-dimensional theory to construct a higher-dimensional action, it is often 
referred to ClS db holographic model. 

To adequately explain the system leading to the correspondence, we must 
understand a few essential ideas of string theory. The fundamental objects of 
string theory are called p-branes, where p specifies the dimensionality of the brane. 
Two important excited 1-branes are the closed and open strings. As the name 
suggests, closed strings form a loop and propagate freely through space, whereas 
open strings have endpoints that must terminate on a objects called Dirichlet p- 
branes, or just Dp-branes. A simple illustration of these objects is shown in Figure 
2.1[ Open and closed strings are characterized by their lengths Is and a coupling 
constant Qs that determines the strength of their interactions. The string length 
is often expressed in terms of the Regge slope parameter a' [3T] . 



Given a stack of parallel D3-branes in a type IIB string theory, two low- 
energy descriptions arise depending on the permutation of the low-energy and 
strong-coupling limits. Taking the low-energy limit of the open string excitations 
first, we see that a SU(A^) gauge theory results, specifically the Af = 4 SYM 
theory with a gauge coupling gvM- The relation between the gauge coupling and 
the string coupling becomes [321 E3] 



Assuming that the string length is small, we neglect the interactions of the closed 




(2.1) 



9ym = ^^9s- 



(2.2) 
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Figure 2.1: A simple illustration of closed strings (loops) and open strings which 
have endpoints on the D-branes. 

strings which exist away from the D-branes in the bulk, that is far from r = oo. By 
taking the strong coupling limit next, we obtain a strongly coupled SYM theory. 

On the other hand, taking the strong coupling limit first means the back- 
reaction from the coincident branes cannot be neglected and the resulting 
curvature of spacetime must be accounted for using general relativity. One finds 
black-hole solutions to the supergravity equations where the classical metric takes 
the form as in [321 [331 El] , 




where the curvature radius R is given by 



(2.4) 



Because of gravitational red-shifting, the low-energy limit corresponds to moving 
away from the D3-branes at the ultraviolet (UV) boundary located aX r = oo. 
Taking the horizon limit when r <^ R, the supergravity metric (2.3) takes the 
form, 



R 



(2.5) 



where r is strictly positive. We can separate the metric (2.5) into S"^, 

R^dVL^, 

and AdS^, 



R 

-dt^ + dx^) + -lydr'^, 



(2.6) 



giving an AdS^ x string theory. We use an alternate form of the AdS^ metric 
where 



_ W_ 
r 

This transformation gives a metric with a single warp factor, 

ds"^ = a{zf [-dt^ + dx^ + dz^) , 



(2.7) 



(2i 



where a{z) = R/z, 2; > 0, with a UV cut-off residing at an infinitesimal value, 
z = Ro- 



As shown in Figure |2.2[ we began with a stack of D3-branes and end with two 
low-energy, strongly coupled systems. In AdS/CFT, classical gravity using the 



metric (2.5) adequately describes the full ten-dimensional string theory as long as 

(2.9) 



the radius of curvature is much larger than the string length 

R > Is. 



Using ([2^, ([2^, and ([2^, we see that 

R^ 2 
Y = 9ymN. 



(2.10) 
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Thus, (2.9) requires that the quantity QymN, often referred to as the 't Hooft 
couphng A, is large, 

(72^iV = A»l. (2.11) 

In other words, classical gravity is valid whenever the corresponding gauge theory 
is strongly coupled. The duality between string theory and the SYM theory should 
not need any corrections from stringy effects. We note here that A plays the role 
of the effective gauge coupling. The gauge coupling is usually tied to a scalar field 
in the string theory called the dilaton, $, where 

$ = logA. (2.12) 



Stack of N DS-branes i > SU(N) Gauge Theory 



strong 
Coupling 
Limit 



Low Energy Limit 



Curved-space 
gravity solution 
for JV D3-branes 
and closed strings 



Strongly 
Coupled Gauge 
Theory 



Closed strings 
in D3-brane 
throat 



Figure 2.2: The order in which the limits are taken cause the stack of D3-branes 
to either look like a strongly coupled gauge theory or a low-energy supergravity 
theory [3]. 



11 

In order to work with this correspondence, we reduce the ten-dimensional 



metric (2.5) to a five-dimensional AdS metric. Eliminating the S manifold on 
the gravity side removes some symmetry on the gauge theory side. The R-charges 
associated with supersymmetry happen to have a S0{6) symmetry, the same 
symmetry as the five-sphere. Removing this symmetry, we effectively reduce the 
SYM theory to a non-supersymmetric CFT. By disregarding the 5-sphere 5*^, we 
reduce the ten-dimensional string theory to a five-dimensional classical gravity 
solution and the A/" = 4 SYM theory to a simpler four-dimensional, conformal 
Yang-Mills (YM) field theory, our AdS/CFT. 



2.2 The Correspondence Dictionary 

While no rigorous mathematical proof establishing the AdS / CFT correspondence 
exists, it passes many nontrivial tests that have lead to the creation of a corre- 
spondence dictionary |35]. In general, the dictionary allows one to draw connec- 
tions between a strongly coupled gauge theory in d dimensions and a gravity dual 
formulated in AdSd+i- There are two methods for constructing the duality: top- 
down and bottom-up. Top-down models start with a critical or non-critical string 
theory and reduce to a lower-dimensional supergravity theory. A correspondence 
is then drawn to a strongly coupled gauge theory. Bottom-up construction be- 
gins with the precepts of a strongly coupled gauge theory and builds an effective 
gravity dual in AdS space. All the work presented in this thesis uses bottom-up 
construction. 

Essentially, the correspondence equates the generating functional of gauge- 
invariant operators of the gauge theory with the minimum of the supergravity 
action, subject to boundary conditions of the supergravity fields coinciding with 
sources [281 [36]. other words, 

(eW«''-</'°(-)c^(-))cFT = e'^suGRAW|^^^^ ^ (2.13) 

where (p^ is a source term residing on the boundary of the SUGRA model. The 
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gauge theory resides on the boundary, dAdS, of the gravity duaL Thus, for every 
operator of conformal dimension A in the gauge theory, we have a couphng term. 
As a result, we insert a scalar field with mass m into the AdSd+i gravity 
dual. Considering the Lagrangian of a free scalar field, 

>C = (?^^9A/05^0-mV, (2.14) 

the momentum space field equation becomes 

z'+% {z^-%(j)) - (m^ + kh^) 0, (2.15) 

where df(j) = k'^4>- Near the z = boundary, the k'^ term is neglected, and the 
field behaves asymptotically as [271 129] 

(f)^(j)^z^- + {0)z^+, (2.16) 

where the two roots are 



A± = -±^-+m2i?2. (2.17) 

In most cases, we specify the operator dimension A and d such that the mass is 
determined by 

m2 = A(A-rf), (2.18) 



giving the two possible solutions of (2.17) 



A+ = ci, A_ = ci-A. (2.19) 

It is also possible, however, to couple an operator with a p-form B, 

BAO. (2.20) 



'Ma 

In this case, the scalar field in the gravity dual must take a mass of [29] 



= {A-p){A-p- d). 



(2.21) 
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At first glance, it seems that we are able to model a limited subset of gauge field 
operators, namely when A > p + d, to avoid instabilities from tachyons. However, 
it has been shown that tachyons above a certain mass, 

> (2.22) 

are allowed [37]. 

In addition, the correspondence dictionary allows for symmetries of the gauge 
field to be mapped to the gravity dual. If the gauge theory possesses a global 
symmetry represented by a transformation t/ = e*^, then rj is a constant on the 
boundary of the gravity dual. In the bulk {z ^ 0), nothing prevents rj from 
becoming a function of spacetime coordinates, x and z. Therefore, the global 
symmetry in the gauge theory becomes a local symmetry in the gravity dual. In 
that case, one inserts a gauge field into the gravity dual for every relevant global 
symmetry. 



2.3 The AdS/QCD Correspondence 

Now that we have addressed the basis of the AdS/CFT correspondence, we will 
discuss how it relates to QCD. One immediately sees problems in connecting the 
CFT to QCD. By definition, the CFT has no intrinsic energy scale and no mas- 
sive particles, whereas QCD has a scale set by Aqcd and a number of massive 
resonances. Similarly, a conformal theory has no running coupling, whereas QCD 
has a coupling which increases its strength as energy decreases, leading to con- 
finement. Before we explain the means of breaking the conformal symmetry and 
incorporating confinement, we summarize some of the characteristics of QCD with 
a large number of colors Nc, the gauge theory that will correspond to the AdS 
gravity dual in Chapters |4] and |5j 
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2.3.1 Large-iVc QCD 

The precise motivation for connecting QCD (or any gauge theory) comes from 
considering the 't Hooft laige-Nc hmit [38j. QCD possesses a SU(3) gauge sym- 
metry, but no expansion parameter. By replacing the gauge group with SU(A'"c) 
and taking the Nc ^ oo hmit, 't Hooft generahzed QCD at infinite colors and 
performed expansions out in powers of l/Nc- The large- A^^c limit is a natural 
extension of a low-energy string theory construction, as we will show 

In large- A^c QCD, in which we assume Nc ^ Nf, gluons dominate the gauge 
theory. When drawing the relevant Feynman diagrams, it is best to consider them 
in double-line notation where one line represents one color charge so that gluons 
are represented by a double line. Each free color index then contributes a factor 
of Nc and each vertex a factor of qym to the amplitude. The diagrams naturally 
fall into two categories: planar and non-planar. Planar diagrams are those drawn 



without crossing lines [39], as shown in Figure 2.3 Each planar diagram scales 
as A*"^ in a power series of A'~^, where / is the number of loops in the diagram. 



Non-planar diagrams, such as the one shown in Figure |2.4[ also expand in powers 
of Nc and A. However, the non-planar diagrams are suppressed by a factor of A^^ 



c 



relative to the planar diagrams in Figure 2.3 



The connection with string theory occurs when associating each Feynman di- 
agram to a Riemann surface. In double-line notation, each line is a closed loop 
defining a boundary of a two-dimensional plane. We construct the Riemann sur- 
face by gluing these planes along their boundaries according to its Feynman dia- 
gram [4j. We find that the amplitude. A, of any given diagram is an expansion in 
terms of N^ and A whose powers are fixed by the genus, g, of the Riemann surface 



A = J2N'c'''m- (2.23) 

9=0 

All the planar diagrams correspond to Riemann surfaces with g=0, whereas non- 



planar diagrams must map to surfaces with g > 0. According to (2.23), diagrams 
associated with surfaces of g > are suppressed by powers of A"^. 
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Figure 2.3: Three examples of the most simple planar diagrams in double-line 
notation j4j. Each free color index contributes a factor of A^^^ and each vertex 
contributes a factor of gvM- 



To summarize, we see that each diagram in the gauge theory is associated with 
a Riemann surface. All non-planar diagrams are suppressed as — )■ oo. Thus, 
the planar limit of the gauge theory corresponds to the classical limit of the string 
theory where only Riemann surfaces of (? = are considered. This means only a 
subset of Feynman diagrams are relevant for the large- iVg limit of QCD [U [221 ES] • 

When considering the large- iVg limit, some characteristics of QCD mesons are 
modified. The decay width of mesons scales as 1/Nc so that a.s Nc oo they 
become infinitely narrow states. It is important to note, however, that for highly 
excited states where n ~ N^, widths can become large. Also in large- A'^c QCD, we 
can express the two-point correlators as 

°° 77-2 

{0{k)0{-k)) = J2^^^ (2.24) 

n=l " 
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Figure 2.4: The most simple non-planar diagram in double-line notation is shown 
[1] . The amplitude of this diagram is suppressed by a factor of relative to the 
planar diagrams. 



where the matrix element for O to create the nth meson of mass m„ is 



(O|0|n). The sum in (2.24) must be infinite since the two-point function itself 
diverges logarithmically for large k"^ [3H1 SO]- The infinite tower of masses cre- 
ated in the large-A^^c limit possesses similar qualitative features of the so-called 
Kaluza-Klein tower. As we show in Chapter [3} Kaluza-Klein decomposition is an 
important method in analyzing the meson mass spectra. 



2.3.2 Conformal Symmetry Breaking and Confinement 

AdS space is conformally invariant; therefore, the string theory corresponds to a 
conformally invariant gauge theory with characteristics of large- A^c QCD. For the 
gravity dual to describe a gauge theory with desired properties, we must break the 
conformal symmetry on the gauge side. Since pure AdS space on the string side 
ensures conformal symmetry on the gauge side, the simplest solution is to modify 
the AdS space in some fashion. It should be noted, however, that QCD becomes 
roughly conformal at high energies, so we should recover AdS space in the limit 
z — )■ 0. Imposing a hard wall or a soft wall are two methods of modifying AdS 
space on which we will concentrate. Breaking conformal symmetry in this manner 
also mimics confinement by restricting the effective range of the gauge theory. 
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The Hard Wall 

The hard wall refers to the simplest means of breaking the conformal symmetry 
of the gauge theory. First investigated in [ITl H2] , further work was performed in 
[SI im SSI SSI Slj- In hard- wall AdS/QCD models, one imposes two branes at 
z = Rq and z = Ri called the UV and IR brane, respectively. The cutoff at Ri 
serves two purposes: (i) introducing a scale, 

Ri ^ j^, (2.25) 

^^QCD 

and (ii) mimicking confinement. 

As a phenomenological model, the hard wall was successful in describing a 
number of QCD properties, such as form factors, effective coupling constants, 
and correlation functions. However, it failed to accurately model meson mass 
spectra. Regge trajectories, the behavior of mass spectra scaling ~ n, is well- 
established experimentally, but AdS/QCD models using the hard- wall geometry 
produced mass spectra exhibiting the scaling rUn ~ n. 



The Soft Wall 

Introduced as a refinement to the hard wall [HI SSI EOl EH |52l [53], the soft-wall 
geometry imposes a z-dependent dilaton field to break conformal symmetry and 
mimic confinement. In string theory, standard decomposition of the closed string 
sector yields a one-particle state called the dilaton. When constructing a effective 
gravity dual in AdSd space, the dilaton enters in the action as 

S = j rf'^xe-2*(^)v/=^/:grav + j rf"a;e-*(^)v^£„ + . . . . (2.26) 

The dilaton modifies the equations of motion; therefore, it modifies meson mass 
spectra, form factors, correlation functions, and any other property derived from 
the action. 

In most works, the dilaton usually has a power solution. 



(2.27) 



IS 

where /i introduces a scale for the gauge theory. For simphcity, /i is often set as 

/i = Aqcd. (2.28) 



Given (2.27), the meson mass spectra behaves as 

ml ~ n^-^ . (2.29) 

If z/ = 2, we recover Regge trajectories in the mass spectra. We comment on this 
case in Chapter |3| 

2.3.3 Running of the Coupling 

At high energies, QCD is weakly-coupled and has a coupling with a logarithmic 
running, 

«s(Q) = /^l^l, (2.30) 

l + laM log^ 

where a<j is proportional to the square of the coupling constant g, and [5l] 

/5 = -yiV,+ ^iV^. (2.31) 

The P function becomes much more complicated at low energies and cannot be 
written in closed form. Experimental evidence, however, has shown that the QCD 
coupling increases at low energies. 

While conformal symmetry is broken in either the soft- or hard-wall model, 
only the soft-wall model allows for the coupling constant to run. In string theory, 
the string coupling often is not constant [31j. The coupling is determined by the 
dilaton field 

^7,iV,~e*(^), (2.32) 

making the 't Hooft coupling A, 

A = e*(^\ (2.33) 
By defining the coupling in this way, we incorporate a running coupling, an ad- 



vantage of the soft-wall model. In general, we could specify (2.33) to mimic the 
running seen in QCD. 
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2.4 Conclusion 

In this chapter, we outhned the precepts of the AdS/CFT correspondence dic- 
tionary, and tied the gauge side of the duahty to a large-A^c QCD theory. Wc 
covered two key means of building a gravity dual model from a gauge theory: (i) 
gauge theory operators map to scalar fields in the gravity action and (ii) global 
symmetries in the gauge theory map to local symmetries in the gravity action. 
In order to deform (break) the AdS space (conformal symmetry) and introduce a 
mass scale, we introduced methods of sharply cutting off the extra dimension with 
a IR brane at a point z — Ri (the hard-wall model) and effectively cutting off z 
using a dilaton background field (the soft- wall model). Running of the couphng 
only occurs in the soft-wall model. 



Chapter 3 

Applications of AdS/QCD 



"History proves abundantly that pure science, undertaken without regard to 
apphcations to human needs, is usually ultimately of direct benefit to mankind." 

-Irving Langmuir 

Now that the basics of AdS/CFT and its phenomenological cousin AdS/QCD 
have been reviewed, we detail some of the models that have been formulated in 
the last decade using the gauge/gravity duality. We comment on models ad- 



dressing hadronic properties in Section 3.1[ We then outline how to incorporate 



thermodynamics into the AdS/QCD picture in Section 3.2 



3.1 Hadronic Properties 

Work addressing phenomenological model^ of hadronic properties are prolific 

[SQlEIlEllHlllalliailHlEIlligilHlllSlI^ 

[63t Ell ESI EHl EZl EH EH]- The models that concern our work usually begin with 
a relatively simple action incorporating a limited number of QCD properties: (i) 
chiral symmetry, (ii) quark and gluon condensates, and (iii) a background dilaton 
field to model a running coupling. 

^ In this context, phenomenological models refers to those that do not dynamically generate 
a background geometry, but assume one. 
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The string-frame action for these models follows the standard AdS/QCD dic- 
tionary. Gauge fields Lm and Rm model the SU(A^/)lX SU(A^/)i? global chiral 
symmetry of QCD with Nf flavors of quarks. We also include a bifundamental 
field X"^ that is dual to the operator q"'q^ and takes on a z-dependent vacuum 
expectation value (VEV) v{z) to break the chiral symmetry. Since X"-^ is dual to 
a dimension-3 operator, its mass becomes 

m\R^ = A(A - d) = 3(3 - 4) = -3. (3.1) 

The field X actually becomes a complex field to incorporate the scalar S and the 
pseudoscalar P fields, 



v{z) 



^ab ^ + SH'j e'^"''. (3.2) 

We only will write the indices a, b of the field X when needed. The five-dimensional 
action takes the form. 



|M|^ + m^W+-i,(F^ + Fl) 



(3.3) 



where is the five-dimensional coupling constant rescaling the field strength 
tensors and Fr, which are defined as 

pMN ^ d^R"" -d^'R^ -i[R^,R% (3.4) 

where = L^'^t"' and Tr[t"t*] = 8°^ 12. The covariant derivative becomes 

D^X = d^X - iL^X + iXR^. (3.5) 

To describe the vector and axial- vector fields, we simply transform the L and R 
gauge fields into the vector (V) and axial-vector (^4) fields where 

= V^ + A^, (3.6) 
R^ = l^^-A^. (3.7) 
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Substituting equations (3.6) and (3.7) into the last term of the action (3.3), 

1, 



Fl + Fl = 2[d^'L^dML^-d''L^dr,L,,+ ^-[L'',L^][LM,L 



N 



2 [ d R OmRn — OR OnRm + - [R ,R ] [Rm, Rn] 



+ 4 ^d'^'A^'dMYN - d^'^A^'dNAM + -[A''\A''][Am, A^] 
= 2(F^ + Fl), 



where now 



T?MN 



QMyN _ QNyM _ ^[V^M^ yN^^ 

a/2 



pMN ^ qM^N _qN^M _^^^M^j^N^ 



(3.8) 

(3.9) 
(3.10) 



Equation (3.8) allows us to rewrite the action (3.3) in terms of the vector and 



axial-vector fields, 

S5 = - j (fx^e-^^'^li 



\DX\^ + ml\X\^ + ^^{F^ + Fl) 



where the covariant derivative now becomes 



(3.11) 



(3.12) 



The 5D coupling constant is inserted in order to match the vector two-point 
function liy^'f) from the operator product expansion of QCD [55] . 

j d''x{j;{x)j'M) = 6''\k,K - gV)ny(A:'), (3.13) 

where k is the 4D momentum and J"' is the vector current. As shown in [Ul 170] , 
the two-point function near the boundary, which corresponds to large scales in 
the field theory, becomes 

1 



2 log/c^ 



(3.14) 
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Matching to the large- A''^ QCD perturbative result, 

^y(A:^) = -^logfc^ (3.15) 

we find that 

d = (3.16) 

To obtain linear mass trajectories from this action, we must introduce a back- 
ground dilaton field with the asymptotic behavior 

0(z ^ oo) ~ /i^z^ (3.17) 

where /i sets the mass scale for the meson spectrum. The ^-dependent dilaton field 



also ensures broken conformal symmetry and a running coupling. In Section 4.2.1 
we show that if chiral symmetry remains broken as excitation number increases, 
indicated by the behavior of v{z), then the quadratic behavior of the dilaton arises 
naturally. 

3.1.1 Equations of Motion 



By varying the action (3.11), we find the equations of motion for the scalar VEV, 
scalar field, pseudoscalar field, vector field, and the axial-vector field. We only 
consider the tree-level terms since we are using a weakly coupled gravity dual. In 
every sector, we derive the mass spectrum in the hard-wall model and the soft- 
wall model. In the hard-wall model, it is the customary to find mass eigenvalues 
from the poles of the two-point correlator ^31 HI]. In AdS/QCD, the two-point 
correlator is found by taking two functional derivatives of the action with respect 
to a source field, ip^, at the boundary. 

In the soft wall, Kaluza-Klein (KK) decomposition is the preferred method for 
finding mass eigenvalues. In KK decomposition, we break the field of interest. 
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ip for instance, into an infinite tower of 4D components and purely z-dependent 
parts, 

oo 

iP{x,z) = J2'^n{t,x)Mz)- (3.19) 



n=0 



Then using Proca's equation [7T] . 

we formulate a mass eigenvalue equation from the equation of motion. 



(3.20) 



Scalar VEV 

Since the scalar VEV only depends on the bulk coordinate z, it does not require 



KK decomposition. We merely use (3.2) and vary the Lagrangian with respect to 
v{z), 



6C 



VEV 



-6 e 



-gTi 



^ A A , 2 ^ 



-6 i^-e-M^f {a{z)-^d,vd,v + m\v') 



finding the equation of motion to be 



a 



div — dzddzV H — —d,v — a^m\v = 0. 



(3.21) 



(3.22) 



For the AdS metric and m^i?^ = —3, we see that the equation of motion reduces 
to 

(3.23) 



" A' ' 3,3 

V — (pv f H — -V = U, 



z z^ 



where (') denotes a derivative with respect to z. 



is 



In the hard wall model where 0' = 0, we find that the exact solution to (3.23) 

Vhw{z) = ciz + C2Z^, (3.24) 
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where Ci and C2 are integration constants. Comparing this solution to the ex- 



pected behavior (2.16), ci and C2 correspond to the source term and the operator 



expectation value, respectively. 



Cl ~ nig, 

C2 ~ (qq) = cr. 



(3.25) 
(3.26) 



However, |l3l HI] argue for an alternative method. Placing a UV (IR) brane at 
z = Rq {z = Ri), they specify the hard wall in the range Rq < z < Ri, where 

R 



rrin 



Ri 



-Vhw{Ro) 



a = RvhwiRi] 



making 



Cl 



C2 



uigRl — aRl 

RR\ {Rf — -^0) 
cr — mqRi 



(3.27) 

(3.28) 
(3.29) 



RR\ {Ri — -^0) 

In the standard soft-wall model where (p = {fJ^zY, the solution to (3.23), is 
given by mi ED] 

VsA^) = '^zT(l)u(lo,f^'zA, (3.30) 



R \2J V2 / 
where U{a, b, y) is the Tricomi confluent hypergeometric function. One possible 
solution is disregarded since the action must be finite in the IR. In the small-^ 



limit, (3.30) expands to [50] 



'Jsw{z) —Z 



(l-27^-21og(^z)-^(^ 



(3.31) 



where is the Euler function. However, comparing (3.31) to the standard form 



(2.16), one sees that the chiral condensate is proportional to the quark mass, in 



contradiction to QCD. This proportionality stems from the fact that only one 



solution survives the imposed constraints. Furthermore, the solution (3.30) has 
an asymptotic limit v{z) — )■ m^/ fi for large values of z, which influence the large-n 
excitations. Since v is responsible for the mass splittings between the p's and ai's. 
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the meson masses become degenerate in the chiral hmit, mg — )■ 0. this does not 
occur in QCD. The problem of couphng between exphcit chiral symmetry and 



A particular normalization scheme for the scalar VEV is implemented in the 
small- 2; limit jSHI, which we use in Chapter |4| The normalization scheme arises 
from the fact that when coupling an operator O to a source 0° one always has 
the freedom to redefine O — t- (^O and (p^ — ?■ (p^/C, keeping the coupling term 0(f)Q 
unchanged [55]. In the case of v{z), we define a new parametrization, 




We fix ( by matching the two-point correlation function of the scalar field between 
the QCD result at large Euclidean momentum, k"^, and that calculated from the 
gravity dual. In QCD, the two-point correlation function is 



spontaneous chiral symmetry breaking is addressed in Chapter |4| 





(3.32) 



/ 



rfW^"(gg(x)gg(0)) 




k^logk"^ + . . . , 



(3.33) 



whereas the dual model yields 



/ 



d^xe''''-'{qqix)qqiO)) 



loge + .... 



(3.34) 



This suggests that the normalization in (3.32) is 



c 




(3.35) 



we assume this relation in Chapter |4j 



Scalar Sector 



Mass eigenvalues arising from the scalar sector represent the /o mesons. The scalar 
field 5* has a similar equation of motion to v{z). The variation of the Lagrangian 
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in terms of S gives 

6Cs = 6{e 



^ Tr [g''''dMS{x, z)dNS{x, z) + m\S{x, zf] ) 



leaving the equation of motion, 



5™2 



a m 



X 



S) 6S, 



e^a-^d, {e-'f'a^d.S) + d^d^'S - a^m\S = 0. 



(3.36) 



(3.37) 



In the hard-wall model, we need to find the poles of the scalar two-point 
function Us- Taking two functional derivatives of the effective action with respect 
to the scalar source field gives the two-point function Ils{k), 



n 



6S, 



(3.38) 



where in this case, 



S^ff = J d^x^g''Sd,S. (3.39) 
After a Fourier transform S{k,z) = J d^xe^^^ S{x,z), it is convenient to define 



making (3.37) 



S{k,Ro) = S\k), 



^d, {a^d,S{k, z)) - k^S{k, z) - a^m\S{k, z) = 0. 



We solve (3.41) directly to find that 

S{k,z) = Ciz^Jiiikz) + Ciz'^Yiiikz), 



(3.40) 



(3.41) 



(3.42) 



where Ji and Yi are Bessel functions of the first kind. In order to obtain a 



zero mode from (3.42), we modify the Neumann boundary conditions to obtain a 



boundary mass term |35j . 



z=Ro,Ri 



(3.43) 
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where ^ is a constant parametrizing the boundary mass. It is convenient to rear- 



range the constants in (3.42) such that 



S{z) = Csz^ (Jiiikz) + hs{k)Yi{ikz)) . 



(3.44) 



Applying the boundary conditions to ( |3.44[ ), we find that 

ikRMikRi) + (1 - j^Ri)Ji{ikRi) 



bs{k) 
Cs 



ikRiYo{ikRi) + (1 - ^Ri)Yi{ikRi] 



RRl{Ji{zkRo) + bs{k)Y,{ikRo))' 
The effective action on the boundary then becomes 



(3.45) 
(3.46) 



Ro 



e 



RR() 



1 + ikRo 



MikRo) + bs{k)Yo{ikRo) 
JiiikRo) + bsik)Y,{ikRo) 



(3.47) 



Taking two functional derivatives of ( 3.47[ ), we are left with the two-point function, 

(3.48) 



TT ^^^ ■UT, M^kRo) + bsik)Yo{tkRo] 



'Ji{tkRo) + bs{k)Yi{tkRo 
In the limit that kRi ^ 1, the poles of II^ occur when 

Jo{ikRi 



Yo{ikRi 



0. 



(3.49) 



The scalar mass m„, where k"^ = — m^, are equal to the poles of the two-point 
correlator. We find an approximate mass tower of 



(3.50) 



where Ri is setting the mass scale in the hard-wall model. We see that indeed the 
hard-wall model fails to produce Regge trajectories. 

In the soft-wall model, where again we take = fi'^z'^, we use a KK decompo- 
sition, 

oo 

S(X, Z) = Sniz)Snix), (3.51) 



n=0 



and Proca's equation (3.20) to simplify (3.36), 
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Sn{z)) Sn{x) 

e'^a-^d, {e-'^a^d,Sn{z)) - {ml + a^ml)Sn{z) 



0, 
0, 



-diS^ + d,<j)d,Sn-3^d,Sn-a'mj,S^-miSn = 0. (3.52) 
In the AdS metric, the scalar mass eigenvalue equation becomes 

~^ 'P ~^ ~o^n — ^n^ny (3.53) 

Z Z'^ 

where (') denotes a derivative w.r.t to z. Following the procedure outlined in 



Appendix |A.2| and substituting 



log 2 



(3.54) 



into (3.53), the equation of motion then simplifies to 

//i/42|3 \ 2 



(3.55) 



Using the form outlined in Appendix A. 1.1, we calculate the scalar eigenfunction 
and mass eigenvalues. 



Sr,\Z 



2^^' ri/ 2 2n rn^ 



2 , 2 



4/x' 



2 ' 



ml = (4n + 6)/x2, n = 0,l,2. 



(3.56) 
(3.57) 



Hence, in the pure soft-wall model, the square of the scalar mass eigenvalues are 
linear in n. We shall compare this result to the results from our modified soft-wall 
model in Chapter |4j 



Pseudoscalar Sector 



The pseudoscalar mass eigenvalues correspond to the pions, the pseudo-Nambu- 
Goldstone bosons. Very few papers actually address the mass eigenvalues of the 
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pseudoscalars HI] . The pseudoscalar sector is the most difficuh to investigate 
because of the many ways to define the pseudoscalar field in terms of the physical 
vr field and because the pseudoscalars couple to a component of the axial-vector 
field. The field P{x,z) has been parametrized as 2tt{x,z) and as tt{x, z)/v{z). 
One can even abandon the exponential representation ( 3.2[ ) and use the linear 
representation [73] . 



In addition, depending on the gauge, the pseudoscalar field couples to or the 
longitudinal component of A^, 

= A^^ + d^^, (3.59) 

generating a system of equations that can only be solved by numerical techniques. 
The various issues surrounding the pseudoscalar sector are addressed in Section 

Vector Sector 

The vector field corresponds to the p mesons. The vector gauge field, Vm is 
introduced into the gravity dual action through the field tensors Fl^r used to 
model the global chiral symmetry of QCD. Equivalently, the vector field is dual 



to the 1-form vector current g7^g. According to (2.21), the mass term rriy. In 
either case, the relevant terms from the action are 



J d:>xe-^^ {dMVNd^^V^ - dMVNd^V^) . (3.60) 



In this extra-dimensional model, we analyze the equations of motion in the so- 
called axial gauge while imposing a Lorentz gauge on the transverse components, 

K = 0, d,vi = 0. (3.61) 
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We drop the subscript ± since we ignore the longitudinal part. Varying the action, 
we find that 



SS 



V 



-S [e-^V^g^'g"'' {d.KdpV, - d.Kd^V,) + e-^V^gg^'g^^'d^V^d^V,) 
-e-'^V^g^'g"'' {d,V,6dpV, - d,V,5d^V,) - e~^,/^g''g^''d,V^6d,V, 

{e-'*'a{z)d,d''V^ + {e-*a{z)d,V^)) (3.62) 



leaving 



(3.63) 



In the hard wall model, just as in the scalar sector, we find the two-point cor- 
relator Ily, whose poles give the mass eigenvalues m„. After a Fourier transform 
V^{k,z) = f d'^x e^''^ Vf^{x, z), we again factor the vector field. 



V,{k,z) = V;!{k)V{k,z), V{k,I^) = l. 



Substitution (3.64) into (3.63), we see that 



a-^d,a{z)d,V{k, z) - k'^V{k, z) = 0, 



whose general solution is expressed as 



V{k, z) = Cvz {Ji{ikz) + hv{k)Y^{ikz)) 



Applying the boundary conditions. 



d,V{k,z) 



z=Ro,Ri 



(3.64) 



(3.65) 



(3.66) 



(3.67) 



we find that 



bv{k) 

Cv 



Jo{ikRi) 
YoitkR,)' 



Ro (JiiikRo) + bvik)Y,{ikRo)) ' 



(3.68) 
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The effective action in the vector sector is 



eff 



d'^xa(z)V^d,V>'\,=Ef> 



Rn'^' JiiikRo) + bv{k)Yi{ikRo) ' 



making the two-point vector correlation function, 

ikR MikRo) + bv{k)Yo{ikRo) 



(3.69) 



Uv{k) 



Ro JiiikRo) + bvik)Yi{ikRo) ' 



(3.70) 



In the hmit as i?o and kRi ^ 1, the poles of (3.70) occur when 

bv{k)Y^{^kRo) = -^|^yi(^fci?o) = 0. (3.71) 
Yo{ikRi) 

The condition looks similar to the one in the scalar sector, and we obtain a similar 
asymptotic behavior for vector mass excitations, 

3\ TT 



m„, ^ |n + - I — , n = 0,l,2, .... 



(3.72) 



In the soft-wall model, we again use a KK decomposition, 

oo 

V,{x,z) = J2yn{z)V,,n{x), 



(3.73) 



n=0 



and Proca's equation (3.20) to find that the equation of motion becomes 



mlVnVn + e^a-^d, {e-^a{z)d,V^) V„ = 0. 
In the AdS metric, the equation of motion simplifies to 

V:-ct>'Vn--Vn + mlVn = Q, 



(3.74) 



(3.75) 



where (') denotes a derivative with respect to z. We transform (3.75) using the 



process described in Appendix |A.2[ where 

T/ _ „, p<f>{z)+\ozz 



(3.76) 
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The equation of motion then becomes 

-f" + ( fi^z^ + ) Vn = mlvn. (3.77) 



4z 



Matching the form in Appendix |A.1.1[ we see that the eigenfunction and mass 
eigenvalues are 



2 



ml = 4(n + l)/i2, n = 0,l,2,..., (3.79) 
where the vector sector shows the same hnear trajectory when plotting ml vs n. 



Axial-vector Sector 

The axial- vector mass eigenvalues should match the mass values of the oi meson. 
The equation of motion for Am is obtained in the same manner as the vector 
field Vm, however, we see the symmetry between the vector and axial- vector is 
broken by the scalar VEV v{z). Now we see why the function v{z) contains the 
parameters for explicit and spontaneous chiral symmetry breaking: it couples to 
the axial-vector field, breaking the symmetry. 

We again use the axial gauge, imposing a Lorentz condition on the transverse 
component, 

A, = 0, 9^v4^ = 0. (3.80) 
Following the same steps as the vector field analysis, we obtain the exact form in 



(3.62) with a symmetry breaking term v A 



SSa = {e-'^a{z)d^A^ + 0, {e-'^a{z)d,A^) + e-'t'glaizfv^A^) 5 A''. (3.81) 
Using the vector sector as a guide, this reduces to a familiar equation of motion, 
e^a-'d, (e-<^a(^)9,A^) + d^d^A^ - g^v^A^ = 0, (3.82) 



where the derivative (') is with respect to z. As f — )■ 0, (3.75) and (3.82) are 



essentially equivalent differential equations; therefore, the vector and axial-vector 
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sector contain a degenerate tower of masses. The scalar VEV controls the mass 
splittings between the two meson mass eigenvalues. 

In the hard-wall model, we again break the axial-vector field up into 



A, = Al{k)A{k,z), A{k,Ro) = l, 



(3.83) 



making the equation of motion, 



.A" + lA' + k^A + gla%^A = 0. 



(3.84) 



The solution to (3.84) depends greatly upon the form of v{z). We are concerned 
with the interesting cases of a constant or linear scalar VEV [IE]; therefore, we 
assume 



vi{z) 



7, 



(3.85) 
(3.86) 



where 7 and F are positive constants. Of course, (3.85) and (3.86) may only be 



the large-2 asymptotic behavior of the VEV, since we have stated that at small 
z, f ~ mqZ + az^. The large- 2; solution corresponds to the asymptotic behavior of 
high excitation modes n since the large-z limit of the effective potential affects the 



higher modes more than the lower modes. The equation of motion (3.84) becomes 



a modified form of the vector equation (3.65). Using fi, the axial-vector solution 
becomes 



A{k, z) = Caz [Ja{ikz) + bA{k)Y^{ikz)] , a = + 



This solution leads to an asymptotic mass spectrum. 



rrir, 



(3.87) 



(3.88) 



Using V2, we find the solution becomes 



A{k,z) = Caz 



Ji (^zi^fk^Vgim^^ +bA{k)Yi (^zi^k'^ + glR^T^ 



(3.89) 



which gives 
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(3.90) 



As expected, the axial- vector mass spectrum is shifted relative to the vector masses 
because of the CSB function v{z). Also anticipated, the hard wall does not produce 
Regge trajectories in the axial- vector sector. A more general treatment of the fir st- 
and second-order corrections to 11^ is done in [72] . 

In the soft wall, we proceed as in the vector sector. We introduce the KK 
decomposition, 

oo 

A^{x, z) = Y, an{z)A^ix). (3.91) 



n=0 



Using the transform in Appendix |A.2| ( |3.82[ ) becomes 



-< + ( + 



(3.92) 



4^2 ^2 

To obtain a solution to and m^, one must again specify the function v{z). 



We consider two forms, (3.85) and (3.86). For the constant VEV (3.85), the 



eigenfunction and eigenvalues are 

an{z) = e-^^' 



2zq\ 



(y/l+gl^^R2 + g)\ 

(^4n + 2^Jl + g|J^R^ + 2^ 



where 



4/i2 



Using a linear VEV (3.86), the eigenfunctions and eigenvalues become 



an{z) 



2z{q + w)\ 1 2 2^ 
(1 + ^ + ^)!%+-^^ 



where 



(4n + 4 + glT'R') /.^ 
glV^R^ 



4/^2 



(3.93) 
(3.94) 

(3.95) 

(3.96) 
(3.97) 

(3.98) 
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One sees that indeed the vector and axial- vector mesons are degenerate if f — )■ 0. 
In addition, their Regge trajectories slopes are parallel, precisely what experimen- 
tal mass values of p's and ai's suggest [1]. 



3.2 Thermodynamics 

Thermodynamic quantities become relevant in the AdS/QCD framework when ad- 
dressing different physical phases of QCD. At a critical temperature hadronic 
matter melts into a strongly coupled plasma where quarks and gluons freely as- 
sociate [m [151 m] • We describe the thermodynamics of this quark-gluon plasma 
(QGP) by modifying the metric of the gravity dual model. The gravity dual met- 



ric (2.8) corresponds to the gauge theory below T^. Above T^, we must consider a 



new metric, 

ds^ = a{zf (^-f{z)dt^ + dx^ + , (3.99) 

with boundary conditions /(O) = 1 and f{zh) = 0. The boundary conditions of / 
allow for a useful parametrization that simplify expressions, 

(3.100) 

where 

r 1 

V{z) = / dx—-^. (3.101) 

Jo a{xr 

The black-hole thermodynamics, outlined in |75l [76|, [77] then play the role of the 
QGP thermodynamics. The position of the black-hole horizon Zh determines the 
temperature, and, ultimately, influences all the thermodynamic quantities of the 
gauge theory. 

We consider energy expressions to derive the black-hole thermodynamics. For a 
normal gas of temperature T, pressure p, entropy S, and volume V, the differential 
change of the internal energy is simply 



dE = TdS - pdV. 



(3.102) 
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Analogously, the energy of a black hole depends on the three quantities an observer 
can measure: charge g, horizon area A, and angular momentum J, 



dE 



8nG, 



-dA + qd{p + Q dJ, 



(3.103) 



where Gd is the gravitational constant in d dimensions, k is the surface gravity, 
the electrostatic potential, and Q the angular velocity [75j. Comparing (3.102) 



and (3.103), we draw a correspondence between the quantities. 



TdS 



-dA, 



-p dV qdif + fl dJ. 



(3.104) 
(3.105) 



For the purposes of our work, we are only concerned with (3.104), where conven- 
tions dictate that 



T 



S 



K 

2' 
A 

AG, 



(3.106) 
(3.107) 



The area of the black-hole horizon is calculated from the induced metric 7, 



A = j dV^, 



(3.108) 



where dV is the induced volume element. Using the metric (3.99), we find that 
the entropy density is given by 

S a(zhY 

where p is the number of spatial dimensions. 

We stated above that all thermodynamic quantities ultimately come from the 



metric (3.99) evaluated at the position of the horizon z^. To see this, we rewrite 



the t- and 2;-components of the black-hole metric in an expanded form around 
as done in 



ds' 



-jt{z - Zh)dt^ + 



Z- Zh 



-dp'^ -f- . . . , 



(3.110) 
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and make the substitutions, 



p^2^^,{z-zn), (3.111) 
t = ir, (3.112) 



making 



dp" = -^^dz''. (3.113) 



In these coordinates, the metric then appears in polar form, 

ds^ = -^pcir^ + + . . . , (3.114) 

with the variable r and p playing the role of the polar and radial coordinate, 
respectively. Examining an arc length, a, in this coordinate system, we see that 
at constant p = Po and a period of /3, 

The arc length should be equal to the circumference C = 27rpo in order to ensure 
regularity at the origin, 

(2npor = -P'^', (3.116) 

1 1 Fh 
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Following the conventions of finite temperature field theory, we equate the period 
in Euclidean time with the inverse of the temperature, 

r=l, (3.118) 

giving us the Hawking temperature expression. 



Using the metric (3.99), the temperature can be expressed as 

dV{z) 



T = - 



Ati dz 



1 
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(3.120) 



AnVizh) dz 
giving a direct relationship between T and Zh- 

Knowing the entropy and the temperature allows us to find other useful quan- 
tities that emerge in Chapter [5| such as the free energy and the speed of sound 
through the QGP-like medium Vg- From thermodynamics, the free energy is de- 
fined as 



SdT. 



(3.121) 



We can express (3.121 ) strictly in terms of the metric. 



J" 



V 



dzh \A7ia(zh)V(zh] 



(3.122) 



The speed of sound through the medium is also dependent upon S and T 

2 _ dP _ dP/dT _ s 
~ ~de ~ de/dT ~ Tds/dT 
d\ogT 



(3.123) 



dlogs 

where P is the pressure and e is the energy density. Equipped with these relations, 
we are able to calculate much about the deconfined phase of the gauge side of the 
gauge/gravity duality. 



3.2.1 Finite Temperature Action 

When considering the thermodynamics of AdS/QCD, we are interested in dynam- 
ically generating the geometry and the fields. Dynamic solutions are not easy to 
find in AdS/QCD models; most work focuses on phenomenological models to de- 
scribe the thermodynamics of the corresponding gauge theory. The gravitational 
action needed to generate the metric and fields has an intrinsic coupling with the 
dilaton, 

S= f d'^a;e-2*v^£g,,,+ f d^xe-^ U^tter- (3.124) 
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The action considered in this section is widely used [SZl ESI EOl [HD El ESI El 
85| [86t [871 [881 \89\ 190] and contains simple elements reduced from a critical string 
theory to five dimensions [82], 



-"grav 



j c/^xe-2*v^ (i? - A + Ag^^du^dN^ - K($)) , (3.125) 



where (3.125) is defined in the string frame. 

We use the gauge/gravity duality to relate the 5D gravitational coupling Gs 
to the number of colors N^. in the 4D gauge theory. First, consider an AdS^ x 
background. The ten-dimensional gravitational constant Gio = Sn^g'^ll, where Qs 
is the string coupling and Is is the string length |82]. After compacting to five 
dimensions we obtain 

G^o = G5V5 = G^TT^R^ (3.126) 
By the AdS/CFT correspondence, we also have the relation p7] . 



It ■ 



(3.127) 



Combining (3.126) and (3.127), we can then relate G5 to the number of colors, Nc 

1 



(3.128) 



167rG'5 87r2i?3- 

This relation is specific to the AdS^ x 5*^ background, but the numerical factors 



can be different in the case of a pure AdS manifold. Analogous to (3.128), we 
define in our five-dimensional case, 



Nl 



(3.129) 



167rG5 

where M5 is the 5D Planck scale and vq is a constant. The constant vq is fixed by 
matching the free energy to that obtained in the high temperature limit of QCD. 
The result is |1] 

{M^f = (3.130) 



457r 



We use (3.130) for our holographic model in Chapter [s] 
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Since we are interested in the gravitational properties of this action, it is most 



useful to put (3.125) in the Einstein frame. Switching from the string to the 



Einstein frame occurs with a conformal transformation, 



givmg 



9mn 



E 



levrGs 

where we make the substitutions, 



g'^''dM(l)dN(l)-VE{(p) 



(3.131) 
(3.132) 

(3.133) 
(3.134) 



We drop the E subscript from this point. In general, the action (3.132 ) contains no 



mass term for the dilaton; therefore, according to the correspondence dictionary, 
the dilaton must be dual to a dimension-4 operator. The simplest dimension-4 
operator is the gluonic operator, Tr(F^,yF'^^). 

In the Einstein frame, we are able to use the unaltered Einstein equations, 

Gmn = SttG^Tmn- (3.135) 
The metric and all the fields of the system are generated from the system of three 



equations, two from (3.135) and one from the direct variation of the action (3.132) 
with respect to the dilaton (p. Using the Einstein metric. 



a[z) 



they become 



-fiz)dr + dx' + 



77 + 3- 



(3.136) 



(3.137) 
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The system becomes extremely hard to solve for all but the simplest cases. We 
examine two cases that have been investigated in [HI]: the conformal case and the 
power solution. 



Conformal Case 



In the conformal case, we set the dilaton to zero. As a result, (3.138) gives pure 
AdS space. 



a{z) 



R 



(3.137) is satisfied when 



f{z) = 1 - 



y4 ' 



and (3.139) holds true for 



V 



12 



(3.140) 
(3.141) 
(3.142) 



This describes a perfectly conformal gas. The temperature is tied to the position 
of the black-hole horizon, 

fizh) 1 



T : 

4:71 TlZh 

Entropy scales with temperature as expected. 



S 



a{Zh) 



^ = N^T 

45ttz? 45 " 



The velocity of sound through this thermal gas is characterized by 
2 _ dhgT 



(3.143) 



(3.144) 



dlogs 
d 



1 45 1 
log . ^ + - log s 



dhgs \3 ^An^N^ 3 



1 

3' 



(3.145) 



The free energy is found from the entropy (|3.144|) 



vr 



SdT = -—N^T\ 
45 ^ 



(3.146) 

Again, we see that a pure AdS gravity dual produces properties of a perfect fluid. 
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Dilaton Power Solution 



Because the system of equations (3.137), (3.138), and (3.139) is difficult to handle 



we define 6 simple power solution. 



(3.147) 



With this dilaton, however, there is no obvious way to satisfy (3. 137), (3. 138) 



and (3.139) with a temperature-independent expression for V{z). A temperature- 
dependent V invalidates any useful traditional thermodynamic relations. Much 
has been written on this point; some writings focus on extremely complicated 
forms for (p and a{z) jHU EHl |91], but a closed-form, temperature-independent 
solution has yet to be found. Currently, only phenomenological models work 
when assuming (3.147). In Chapter [sj we attempt a novel approach to generating 
a dynamic solution in the soft-wall model. 



Given (3.147), we have 



(3.148) 



which gives a temperature of 



At^zk 1 - r 



— tl^ 7^ 



(3.149) 



We see that (3.149) has a minimum and has two distinct asymptotic limits 

1 



T{ZK -> 0) 
T{zh -> oo) 



T^Zh 



^ .-1 

At, ^ ■ 



(3.150) 
(3.151) 



Three distinct temperature-horizon relationships arise depending upon on the 



value of V as shown in Figure 3.1 The v < \ case is effectively equivalent to 
the conformal case, which was already discussed. The v = \ case has the T{zh) 
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approaching a constant as ^ oo. Confining solutions exist for z/ > 1 [HI]. Tliey 
contain a minimum temperature below which no black hole exists. Above the 
minimum temperature a stable black hole and an unstable black hole can form 
depending upon the slope of the T vs Zh plot, 

stable, (3.152) 

unstable. (3.153) 



dT 

dzh 
dT 

dzh 




Figure 3.1: Depending on the value u takes, the temperature and the horizon 
position have 3 distinct relationships. The v < 1 case produces conformal results. 
When z/ = 1, we find that temperature approaches a constant as, z^ oo. When 
z/ > 1, we find a stable and unstable black-hole solution. 

Stability issues become obvious as we examine the entropy, speed of sound, 
and the free energy for confining theories describing the QGP-like fluid. The 
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Bekenstein entropy is 



3 



1/ ^nu 



4iV2 

S=—^- 3 . (3.154) 

457r 



The entropy behaves quite differently depending on the region of 



^(z.^O,T) ^ -^X^'-^^'-')^ (3-155) 

5(^,^oo,T) ^ , (3.156) 

457r (rT)^ 

where r = Ati/v^'^. As — )■ 0, the entropy goes to the conformal behavior; 



however, the entropy of the unstable black hole (3.156) exponentially decreases 
with rising temperatures. From the temperature and entropy expression, the 
velocity of sound becomes 

S{zn)T\zn) 



S'{zh)T{zh) 



(6 + Qu^,'^z'0 (1 - r (i, -^,^^^zr)) 6 + 6z.;.^< 

For the stable black hole, we see that Vg reduces to 



(3.157) 



.J(.,.-0) « (3.158) 



3 ^ 6 ItttJ 



(3.159) 



For the unstable black hole. 



AttT 1 



jyjJ' a/6' 



(3.161) 



the speed of sound grows without bound as the temperature increases, an unphys- 
ical result. From the Bekenstein entropy (3.155) and (3.156), we calculate the free 
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energy as 

^(..^O.T) = + . (3.162) 

^(..^oo.T) ^ -^|l^l-lir(^i^,(.r)A). (3.163) 

The stable black hole again shows the conformal limit as T increases, and the 
unstable black hole has vanishing free energy with increasing temperature. It is 
easily seen that the thermodynamics reduce to the conformal limit as ^ 
and T — 7- oo, and the unstable black-hole results appear to be of little use for 
describing gauge-field thermodynamics. 



3.2.2 Phase Transition 

The analysis of the gravity dual actions yields information about the Hawking- 
Page transition [921193]. When two solutions exist for a given action, this transition 
occurs as one solution becomes energetically favored over the other. According to 
the correspondence dictionary, the Hawking-Page transition marks confinement 
temperature, Tc, where the phase transition from hadronic matter to a QGP state 
occurs. To calculate the transition, one must consider the on-shell action using 



the metric (2.8) corresponding to a confined state below (thermal AdS), and 



then the metric (3.99) corresponding to the state above (black-hole AdS). A 
simple calculation of the phase transition in an AdS/QCD model was investigate 
in [80], using the free energy as defined in |8 4j , 

= lim [Sbh{Ro) - SthiRo)] , (3.164) 

i?0-S>0 

where again T = 1//3. The transition occurs when J-" = 0. We review the con- 
finement temperature calculation done in a hard- and soft-wall phenomenological 
model [80]. 

In the hard-wall model, we analyze the action 
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where the curvature of either the thermal or black-hole metric is 7^ = — 

In thermal AdS, the on-shell action can be expressed as 

1 r r-^zh rRi op3 

' d^x dr (3.166) 



^hw,th 



where Rq (Ri) is the UV (IR) brane position and J d?x = V. In order to match 
the two metrics at the boundary of AdS, we redefine dt — )■ a/ f{Ro)dt. The action 
becomes 



R'V Rl \ ( \ 



' SG^^^ \C 2zlJ \R\ R^ 
R^V f I 1 Ri 1 



Switching to the black-hole action, we must now adjust the limits of the dt inte- 
gration. The black hole no longer terminates the space at Rq but at Zh, so that 
the black-hole action becomes 

R^V r^'' T'^ 1 
<Shw,bh = z^-^ I dr I dz 



8G, 







dr 






'ro 


1 


1 




R^ 



2'kG^ Jn ./»„ Z^ 

R'V 



Zk(-j-^ . (3.168) 



Using (3.164), we see that the free energy is expressed as 

RsV /J 1 

SttGs \2zf^ R\ 



^ = -^\—,-m]^ (3-169) 



where the transition occurs when 2zf^ = R\. Using the hard-wall relation between 
temperature and horizon position (3.143), we find the critical temperature to be 

m 

Tc = -TT- (3-170) 



In the soft-wall model, we again follow the analysis done in [80] and work with 
a string frame action where $ = ^^z"^, 
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As in the hard wall, the temperature at which the free energy of the system 
vanishes marks the transition point, 



lim {Ssw,bh{Ro) — <Ssw,th{Ro)) 
Ro^O 



lim 

i?0-5-0 



d X 



dt 



dz e — - — 



Ro 



z^ 



d^x / dt a//(-Ro) / dze 
Jo J Ro 



lim 



R^V 
R^'V 

-Zh 



T^Zh 



2G; 



dz — — 

i?.o ^ 



-8R' 



dz 



Ro ^ 
oo 



1 - 



-^0 



dz 



JRo 



dz 



= lim 

i?o-s>0 

R^V 
8G,zl 

where Ei(x) is the exponential integral function. The free energy is then 



e-'^'<{^^'zl-l)+^,'ztEi{-^,'zl) + 



(3.172) 



R^ttV 



e-^ H^^'zl-l) + E\{-^^'zi) + 



8G,zl 

The Hawking-Page transition occurs when 



e-^(x-l) + - + x^Ei(x) = 
2 



(3.173) 



(3.174) 



is satisfied. Since (3.174) cannot be solved analytically, we have defined a new 
variable, x = fi'^z'^, to ease the graphical analysis. We find that the transition 
temperature is [80] 

Te = 0.491728/i. (3.175) 

Again, we see that quantities in the soft-wall model scale with the parameter /i. 
whereas the inverse of the IR brane position I/-R0 sets the scale in the hard- wall. 

The above transition temperatures found above are purely phenomenological. 
The solution used in [HO] was not dynamically produced. Considering the soft- 
wall model in the complete manner, one would need to add a boundary term 
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to the on-shell action called the Gibbons-Hawking (GH) term. It was shown in 
[83] that including the GH term in the calculation above produces a transition 
temperature at T = 0, describing no confined phase. In Chapter [5| we consider 
a dynamic soft-wall solution in an attempt to rigorously calculate a transition 
temperature. 

3.3 Conclusion 

In this chapter, we reviewed the concepts needed to understand the original work 
in Chapters |4] and [5j We summarized much of the relevant literature involving 
meson mass spectra in the scalar, pseudoscalar, vector, and axial- vector sectors 
in both the hard-wall and soft-wall AdS/QCD. Next, we outlined the thermody- 
namics of AdS/QCD used to describe the deconfined state of QCD matter, finding 
the temperature, entropy, speed of sound, and free energy of the thermal fiuid. 
We concluded with showing how to calculate the transition temperature from the 
on-shell actions. 



Chapter 4 



Zero Temperature Dynamics of 
AdS/QCD 



"Sometimes the questions are complicated and the answers are simple." 

Dr. Seuss 

4.1 Introduction 

This chapter works through the dynamics of an AdS/QCD model at zero tem- 
perature. We modify the existing soft-wall version of the AdS/QCD model |1H] 
in order to incorporate several phenomenological features of QCD. First, we in- 
corporate separate sources for spontaneous and explicit chiral symmetry breaking 
(CSB). Second, we show this model does not restore chiral symmetry in highly 
excited resonances [30l . Finally, we produce a new dilaton background profile 
and retain the linear mass trajectories in the scalar, pseudoscalar, vector, and 
axial-vector sectors that make this phenomenological model attractive and fertile 
ground for further research. 



This chapter is separated into five sections. In section 4^, we introduce the 
modified soft-wall Lagrangian. We describe how the new dilaton profile and 
higher-order terms lead to a model describing spontaneous and explicit CSB. In 
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Section 4.3, we explore the QCD meson mass spectrum in the scalar, pseudoscalar, 
vector, and axial-vector sectors of the model. In Section [44 we determine other 
dynamics derived from this theory, such as the pion decay constant, the Gell- 
Mann-Oakes-Renner relation, the coupling of the vector mesons to the pions, and 
the pion electromagnetic form factor. We conclude with a discussion in Section 

m 



4.2 The Dual Model 



We consider a modified version of the soft-wall AdS/QCD model first introduced 
in HH] and further investigated in [21 El Ell EQl [551 |52l [531 [Ml ill ESI [9^ In the 
string frame, the background geometry is not dynamically generated but assumed 
to be five-dimensional AdS space with a metric 



R/z is the warping factor and the Minkowski metric. 



(4.1) 



(4.2) 



where a{z) ^ 

/ -1 \ 
10 
10 
\ 1 / 

The coordinate z has a range < 2; < oo. The action follows from the standard 
AdS/QCD dictionary. We begin with the basic action incorporating CSB with a 



scalar field X dual to the dimension-3 qq operator introduced in Section 3.1 The 



primary shortcoming of the action (3.11 ), however, is the explicit and spontaneous 



CSB entanglement. No chiral limit exists in this theory. As suggested by [38], we 
add a quartic term to the action, whose strength is controlled by the parameter k. 
In section [4.3.2 we show that X actually becomes a complex field to incorporate 
the pseudoscalar field vr. The five- dimensional action takes the form, 

1 



d^x 



DX ' + mt X ' - k X 



2gl 



, (4.3) 
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where k is a constant that will be determined and g1 = /Nf. as shown in 
Section 13. 1[ 

To obtain linear mass trajectories from this action, we assume a soft-wall 
background dilaton, 0, with the asymptotic behavior 

(t){z oo) ~ /i^z^, (4.4) 

where sets the mass scale for the meson spectrum. The 2;-dependent dilaton field 



also ensures that conformal symmetry is broken. In Section 4.2.1 we show that 
assuming that if chiral symmetry remains broken as excitation number increases, 
as evident through the behavior of v{z), then the quadratic behavior of the dilaton 
arises naturally. 

4.2.1 Bulk Scalar VEV Solution 

The field X, which is dual to the operator gg, obtains a 2;-dependent vacuum 
expectation value (VEV), 

Specifying the elements of the bifundamental field fixes the number of flavors 



that we consider. In the case of (4.5), we have two flavors of degenerate quarks. 



The VEV breaks the chiral symmetry SU(2)2,x SU(2)r — SU(2)y by coupling 



to the axial- vector field. Assuming (4.5), we obtain a nonlinear equation for the 



VEV v{z) using the variation principle 

5S = -2e-^^Tr{g^''dMXdM5X + {A,X}{A,5X} + [V,X][V,5X] 

+m\X5X - 2kX^X\X\5x\ = 0. (4.6) 



Taking the trace, integrating by parts, and using |X| = v{z) + S{x, z), we find the 
a simplified expression for the variation. 



6S 
6X 



X^v(z) 2 



^ (d, {e-^^gg''d,v) + e"^ m\v + e'^^g'^v'') . (4.7) 
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Using the metric (4.1), we find the equation of motion for v{z), 

,v{z)) - a^e 



which simphfies to 



K 



v" -[(j)' + -\ v 



mxv[z) - -v[z) 



m\v{z) + —v^ = 0, 



0, 



(4. 



(4.9) 



where (') denotes a derivative with respect to z. 

As noted in [01], highly excited mesons exhibit parallel Regge trajectories 
indicating that chiral symmetry is not restored with increasing n. In order to 
incorporate this behavior the scalar VEV, v{z) must behave linearly as z becomes 
large, 

v{z oo) ~ z, (4-10) 
causing the mass difference between vector and axial-vector resonances to ap- 



proach a constant as z — > oo as shown in Section |3.1.1[ By including a quartic 
term and requiring linear asymptotic behavior for v, we aim to incorporate these 
QCD-like characteristics into the soft-wall model. 



Instead of solving for v{z) directly from (4.9), we assume the VEV asymptot- 
ically behaves as expected, namely 

CR 

(4.12) 



v{z -> 0) 
v(z — > oo) 



^gC , - 3 
Tz. 



Solving for the dilaton (f){z), we see 

1 



<P'{z) 



a^v'{z) 



d,{a%'{z)) - a%mj,v{z) - -v'{z)) 



(4.13) 



where the prime (') denotes the derivative with respect to z. Given the required 



behavior (4.11) and (4.12), we can uniquely determine the dilaton profile up to a 



constant. With this procedure the two sources of chiral symmetry breaking decou- 
ple while simultaneously allowing for linear trajectories in the meson spectrum. 
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A particularly simple parametrized form for that satisfies (4.11) and (4.12) 



IS 



v{z) = -{A + BtanhCz'^ 

1j 



(4.14) 



where A, B, and C are all positive coefficients dependent upon rriq, a, N^., and k. 



Expanding (4.14) at small and large z leads to the asymptotic forms 



v{z^O)L = Az + BCz^ + 0{z^), 
v{z ^ oo)L = {A + B)z. 



(4.15) 
(4.16) 



When A = 0, corresponding to a zero quark mass, the coefficients of the cubic 



term in (4.15) and of the linear term in (4.16) are nonzero, implying a nonzero 



chiral condensate and non-restoration of chiral symmetry. Alternatively, when 
-B = (or C = 0), corresponding to a zero chiral condensate, the coefficients of 
the linear terms in ( |4.15 ) and (4.16) are both nonzero, implying a nonzero quark 
mass and non restoration of chiral symmetry. Thus the parametrized form in 



(4.14) allows the sources of spontaneous and explicit chiral symmetry breaking 



to remain independent. Substituting (4.14) into (4.13) leads to the following 
asymptotic behavior for the dilaton. 



0(;2^O) = -A^z^ + 0{z^), (4.17) 
0(z^oo) = "^{A + BYz^, (4.18) 

where we have chosen the boundary condition 0(0) = to recover a pure AdS 
metric in the z — )■ limit. To reproduce the limits (4.15) and (4.16) the dilaton 



profile at small z (4.17) must differ from that at large z (4.18). This does not 



sacrifice the linear trajectories which (as will be shown) depend on the dilaton 



having the asymptotic form (4.4). Note that the quartic term with strength k is 



necessary to obtain the required behavior. Therefore, modifying the dilaton and 
introducing quartic interaction terms in the Lagrangian is necessary to improve 
the soft-wall version of the AdS/QCD model. 
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From (3.35), the normalization parameter, (, is set conveniently to 



C-P,-^, (4.19) 

(2vr) g5 

where for this phenomenological model, we take Nc = 3. Then the parameters 7, 
A, B and C can be expressed in terms of the input parameters m^, a, /i, 



4 u? 

^ (4.20) 



A = (4.21) 

B = T-A, (4.22) 

C = (4.23) 
VSB 

The input parameters are determined as follows. In order to find /i, a straight 
line is fitted to the versus n plot for n > 3 for the scalar, vector and axial- 
vector mesons, assuming the same slope 4/i^ but different intercepts, namely = 
Afi^n + ttIq. The fit was limited to the higher excitation modes because Regge 



behavior is not observed in the lower modes. Figure 4.1 shows the plot, which 
gives, in units of GeV^, 

/i^ = 0.1831 ±0.0059, 
m|o = -0.6634 ± 0.0038, 
mlo = 0.0806 ± 0.0104, 

m^ o = 1.5023 ± 0.0366. (4.24) 

The quark mass, quark condensate, pion decay constant, and pion mass are all 
related by the Gell-Mann-Oakes-Renner relation, /^m^ = 2mq(j. This relation 
holds in this model as a natural consequence of chiral symmetry [45|; see also 
Section [441 We use the measured values of /tt = 92.4 MeV and m.„ = 139.6 MeV, 
and adjust the quark mass to reproduce the input value of f-,, from a solution to the 



axial-vector field equation in Section 4.4 for a given value of n. The parameter k 



essentially controls the mass splitting between the vector and axial- vector mesons. 
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It is determined to be k = 15 by a best fit to the radial spectra of the axial-vector 
mesons. This results in rUg = 9.75 MeV and S = (204.5 MeV)^. The inferred 
value of the quark mass is consistent with an average of the up and down quark 
masses as summarized in the Review of Particle Physics [T| as appropriate at the 
hadronic energy scale. 




Figure 4.1: A straight-line fit to the measured scalar, vector and axial-vector mass 
spectra for n > 3 used to determine the dilaton mass parameter /x. 



Other parameterizations for v{z) that lead to qualitatively similar behavior 
and asymptotic limits (4.11) and (4.12) exist. Among many possible alternatives, 



v{z) could take the forms of a Pade approximant or a Gaussian, 
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Vp = {ai + a2Z^) / {1 + a^z"^) , 
vg = biZ -b2zexp{-b3z'^). 



(4.25) 
(4.26) 



We explore the alternative parametrizations in Appendix |Bj These forms were 



all studied but the best results were found using (4.14). The parameterizations 



(4.14), (4.25), and (equgaussinO) are compared in Figure 4.2, The resulting plots 



of d(f)/dz and (f){z) are shown in Figures 4.3 and 4.4 It becomes quite evident from 



the figures that a small change in v{z) parametrization leads to drastic change in 
the behavior of the dilaton (t>(z). 



4.3 Meson Mass Spectra 

Our soft wall model can be used to fit the meson mass spectra. It is interesting 
to see how well this simple five-dimensional model matches real data. The scalar, 
pseudoscalar, vector, and axial-vector mass resonances used in our fits are given 
in Tables 4.1[ 4.2, 4.3, and 4.4, respectively. Before proceeding, we give a short 



justification for including selective radial excitations in our analysis. All but one 
resonance of the included states are listed in the Review of Particle Physics [1]. 

The scalar mesons are the most ambiguous particles in this analysis. Un- 
certainty lies in the scalar mesons since mixing is expected among light quark 
mesons, four quark states, ss mesons, and glueballs. This could shift the masses 
of the "pure" radial excitations of the lightest scalar meson by (9(100 MeV). 
There is evidence suggesting that /o(600) and /o(980) contain 4-quark states and 
that /o(1370), /o(1500), /o(1710), and /o(1770/1790) (most recently referred to as 
/o(1790)) are all superpositions of qq and glueballs [97l|98j. However, the 2-quark 
character of /o(980) appears dominant [99j. In addition, /o(1500) and /o(1710) 
most likely have significant glueball contributions because of their absent 77 cou- 
pling, narrow width, enhanced production at low transverse momentum, and the 
relative close proximity to the state /o(1370) and /o(1790), respectively 1100] . 
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Figure 4.2: A plot of v{z)/z for various parameterizations fitted to the mass 
spectra. The best fit to the mass spectra is obtained with (4.14). 



Experiments suggest that /o(1790) makes a natural radial excitation of /o(1370). 
Almost nothing is said of the strength of the case for the /o(2020), /o(2100), 
/o(2200), or the /o(2330). However, /o(2103 ± 8) and /o(2314 ± 25) in [T] agree 
well with /o(2102 ± 13) and /o(2337± 14) [97J, respectively. Therefore, /o(2100) 
and /o(2330) complete our choices for the scalar mesons. 

The pseudoscalar mesons do not have as many intricacies as the scalars. The 
ground-state pion has been known for decades. The excitation modes of the vr 
have been well established as 7r(1300) and 7r(1800) [Ij. We will also include some 
unconfirmed mass states for further comparison. 
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Figure 4.3: A plot of (f)'{z)/z derived from the various parameterizations of v{z). 



The best fit to the mass spectra is obtained with the tanh form (4.14) 



Three vector mesons are confirmed and require httle debate: p(775), p(1450) 
and p(1700) [I]. However, it has long been known that the p(1465) is too massive 
to be the first radial excitation of the p(775) [97J. Studies of the process p + p — )■ 
27r+ + 27r~ suggests the first radial excitation of the p to be 1282±37 MeV, which 
is the value used in our fits. There is evidence of a p(1570), but this is believed to 
be a OZI-suppressed decay mode of p(1700). Based on the rating system in the 
review [97], we include p(1900) as a confirmed mass resonances, and p(2150) and 
p(2265) are included as unconfirmed resonances. 
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Figure 4.4: The dilaton profile (j){z) resulting from the various parameterizations 
of v{z). The best fit to the meson spectra occurs with the tanh parameterization 
(4.14). For z < 1 the behavior deviates from the quadratic asymptotic form (4.4). 



Luckily, the axial-vector mesons are limited in number. While the only con- 
firmed axial- vector resonance is ai(1260), there exists a clear radial progression: 
ai(1640), ai(1930), and ai(2095) [Ij. Two unconfirmed resonances continue this 
progression, the ai(2270) and ai(2340), which are rather close together assum- 
ing Regge trajectories. However, it is clear that ai(2270) is backed with more 
experimental evidence [1]. 
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4.3.1 Scalar sector 



Introducing a quartic term in the Lagrangian causes the scalar excitations to 
couple with their own VEV, giving a modified equation of motion unlike those in 
j. Assuming X{x,z) = {v{z)/2 + S{x, z))e^^ , where P will be defined in 



Section [4.3.2[ and a KK decomposition, 

oo 

S{X,Z) = ^Snix)Sn{z) 



(4.27) 



n=0 



we obtain by variation of the action (4.3), 



-ah-'f'ml^S^iz), 



(4.28) 



where Sn{z) are the Kaluza- Klein modes. Note that by ignoring the nonlinear 



terms in (4.28) we are assuming infinitesimally small amplitudes S'„. Because 



of the z-dependent mass term, (4.28) is difficult to solve analytically for the 



parametrized solution of v{z); however, we implement a shooting method dis- 
cussed in Appendix |A. 3.1 



The scalar equation of motion (4.28) can be brought into a Schrodinger-like 



form using the process in Appendix A. 2 with the substitution 



S„\Z 



Us = (p{z) + 31og^;. 



The eigenvalue equation becomes 
1 



-00 



12 



^11 "J -'^ 2/ \ 



Sn{z) 



where the boundary conditions are 



lim Sn{zo) = 0, 

20-s>0 



dzsJz oo) 



0. 



(4.29) 



(4.30) 



(4.31) 



We find the scalar mass spectra listed in Table 4.1, and displayed in Figure 4.5 



The reproduction of the experimentally measured masses is reasonable, apart 
from an overall constant. This could be a failure of the model in describing this 



62 

sector. However, considering that these hght quark/antiquark radial excitations 
mix with scalar ss excitations, scalar glueballs, and possible four quark states, it 
may be that either the lowest or first radially excited state has been misidentified. 
Removing either the /o(550) or the /o(980) would shift all the higher masses to 
the left by one unit of n, resulting in a much better fit to the model. An obvious 
extension of this work would be to include strange quarks and glueballs and to 
determine the mixing among the resulting scalar states. 




Figure 4.5: Comparison of the predicted scalar mass eigenvalues using the tanh 
form (4.14) of v{z) (solid) with the QCD /o scalar mass spectrum 
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n 


/o experimental (MeV) 


/o model (MeV) 


1 




799 


2 


980 ± 10 


1184 


3 


1350 ± 150 


1466 


4 


1505 ± 6 


1699 


5 


1724 ± 7 


1903 


6 


1992 ± 16 


2087 


7 


2103 ±8 


2257 


8 


2314 ± 25 


2414 



Table 4.1: The experimental and predicted values of the scalar meson masses. 
4.3.2 Pseudoscalar Sector 

Analyzing the pseudoscalar sector is particularly intriguing because of the pseu- 
doscalar representation used in [SI 03] appears to be incompatible with the rep- 
resentations used in jl5l |62l O [60]. Both are in the exponential form 

X = Se'^, (4.32) 

where P is the pseudoscalar part of a complex field X. The former papers de- 
fine P using the vacuum expectation value v{z) of the scalar field, S{x,z), as 
P = 'n-"'{x, z)t°'/v{z), while the latter papers specify P = 271"" {x, z)t°'. The two 
representations produce seemingly different equations of motion and potentially 
different mass eigenvalues. Alternatively, one may use a linear representation, 
X = Xq + i'K^'t"', as in |[73j|. Any physical observables calculated from this model 
should, of course, be independent of the choice of representation. 

In this section, we derive and explore the equations of motion of the two rep- 
resentations, showing that the sets of equations are, indeed, equivalent. One can 
simplify the equations of motion in two different ways: (i) by eliminating one field 
as in [2] , or (ii) by transforming the full system of second-order differential equa- 
tions into Schrodinger-like form. Numerically calculating the mass eigenvalues, we 
find that the eigenvalues found from one method do not match those of the other 
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in our case. Apparently, [2] has found the eigenvalues not of the pseudoscalars but 
of a derivative field, and they are not necessarily equivalent. However, it should 
be noted that the original authors of |2] again assert that methods (i) and (ii) give 
equivalent eigenvalues |101] : we cannot confirm this. 

Concentrating on the pseudoscalar sector, we explore two representations for 
the field X, 

Xe = (t;(2)/2 + ^(x,2))ie2'""("'")*\ (4.33) 
Xi = {v{z)/2 + S{x,z)) l + nr1{x,z)t'', (4.34) 

where 1 is the 2x2 identity matrix. We refer to X^ as the exponential representa- 
tion and Xi as the linear representation. The exponential representation is used in 
[21 El HSl Eni |62] , where it is assumed to be the canonically normalized pion field, 
TTe = Ti/ fn- The linear representation has been specified before in [73], where the 
TT field carries the same dimensions as other fields in the Lagrangian. We already 



see an apparent difference between representations; (4.34) allows for an explicit 



quartic term in vr when substituted into the action (4.3), whereas there is no such 



term in the case of (4.33). The consequence of such quartic terms in vr will not be 
addressed here; however, the quartic term strength impacts the pseudoscalar mass 
spectrum through the parameter k. We only consider field terms up to quadratic 
order. 

Representations 

In this section, we derive the equations of motion arising from the two represen- 



tations (4.33) and (4.34). The pseudoscalar and longitudinal components of the 
axial-vector field mix in the Lagrangian; therefore, we find two coupled differential 
equations for each representation. This makes the numerical work more involved 
than for the scalar, vector, and axial-vector sections, which were already stud- 
ied in [20] • The last part of this section shows that the two sets of differential 



equations are equivalent. 
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Exponential Representation Let us take (4.33) and substitute it into (4.3) 



where we focus only on the terms involving the field 7r{x, z), 

gMPgNR 

H ^ — {dMANdpAn - duANdRAp) 



9l 



(4.35) 



We work in the axial gauge, A^ = 0, and define A^ = Af^± + d^(p, where d^A'^ = 0. 



Separating (4.35) explicitly into four-dimensional components and bulk terms, we 
obtain 



1 



(^) 



{df^iApdyA^ — df^iApdfjAy) 

95 



+ 



^9''9^'' 
9l 



{d^Apd^A, 



(4.36) 



Keeping only terms of the longitudinal part of A^ gives 
1 



-g g^'^iv^df^Trd^'K + v'^d^^d^^ - 2v'^d^-ndy^) 



+V^g''v'd,TTd,TT + 



-gg^'^g^'^ 
~9l 



(4.37) 



Varying (4.37) with respect to tt gives 



5Ce = d,e-^y/^g''v'^d,T:5Ti + e-'''y/^v'^g>'''d^d^{-K - y?)57r. 
Using a Kaluza-Klein decomposition, 

■k{x,z) = ^n„(a;)7r„(z), 

n 
n 

and 



(4.38) 
(4.39) 

(4.40) 
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we can express the system of equations in terms of its 2;- dependent parts. We 
obtain the first equation of motion, 

e<^ 



d, (e ^v'^a^d^'Kn) + (7r„ - = 0. 



(4.41) 



Varying (4.37) with respect to (p and breaking it down into KK modes gives the 



second equation of motion, 



dz'-Pn + (TTn - V'n) 



(4.42) 



Alternatively, we can express (4.41) and (4.42) in a Schrodinger-like form 



We rewrite them in the same form as (A. 15) and (A. 16) in the Appendix by 
substituting 



TTr, 



/(^) = 0(^)+lo, 
9{z) = 



z) + log z, 



(4.43) 
(4.44) 



which eliminate terms involving the first derivative of the fields, vr^ and (p'^. After 
reverting back to the notation 7C it and if) ^ (f the equations of motion become 

0'2 0" (t)'v' 30' 15 



3v' v" 
— + — 

VZ V 



Z^ 

k'2 All 



3 gyn^ 



2z 4^2 



(4.45) 
(4.46) 



where (') indicates the derivative with respect to z. 



Linear Representation When considering the linear representation of the 



pseudoscalar field (4.34), we find quadratic and quartic vr terms that were not 



explicitly present in the exponential representation. After making the appropri- 
ate substitutions, we find that 
1 



2 2 

t> TT 



+ g^''v^d^^d,^ + 



9l 



-dzdf,(pdzdu(p 



(4.47) 
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Once again, we derive two coupled equations. Varying with respect to (j) produces 
a result similar to Xe with the exception of factors of the VEV in the mixing term, 
giving 

e<^5, {^^d.^^ + (vTn - v^r.) = 0. (4.48) 

Varying with respect to vr gives the second equation of the linear representation. 



m 



2^2 



X 



Tin + mlTln = mlvipn- 



(4.49) 



We can express (|4.48|) and (|4.49|) in a Schrodinger-like form with the substitutions, 

(4.50) 

iPn = &ipn, g = - + -iog—. (4.51) 



1, z 

q = 1 — log — . 

^ 2 2 ^ R 



Simplifying the equations and reverting back to the notation tt^ — )■ 7r„ and (pn — )■ 
(Pn for simplicity, we find 

giRv 



-7r„ + 





0" 


3 


0' 

+ 2^ + 






2 


^ 4^ 








3 


30' 


Ki?2^;2 


4 


2 


^ 4^ 


+ 2^- 


2;z2 



(4.52) 



- m„ ) 7r„ = -m„— V3„.(4.53) 



Representation Equivalence The pseudoscalar field representation should 
not affect the physical results obtained from the model. Examining the two sets 



of coupled equations in each representation, we see that neither (4.41) nor (4.42) 



contains an explicit dependence on k, whereas (4.49) does. Although k does not 
appear explicitly in the exponential representation, the dependence enters through 
the function f (-z). 

We begin by expanding Xg, 



X, 



S\ (l + 2i7re + ...) 



+ S* + iVTef . 



(4.54) 
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Comparing (4.54) to (4.34), we surmise that vTef (2;) — )■ tti relates the two repre- 



sentations. Let us substitute vTe — )• 7ii/v{z) into the equations of motion of the 
exponential representation and attempt to obtain the equations of motion of the 



linear representation. The substitution into (4.42) is trivial; it yields 



z z-^ 



(4.55) 



which is equivalent to (4.48) as expected. Showing the equivalence of the other 



two equations requires a bit more work. First, we substitute for TTg in (4.41) and 
then simplify the expression, 



+ ml{ni - vip) = 0, 



(4.56) 



which becomes 



TT, 



' v" — (p'v' v' ) + m^{TXi — Vip) = 0. 



TT, 



(4.57) 



Substituting (4.9), which does not depend on the representation, into (4.57), we 
find 



vr, 



- + 1 vr, + ( - + ) vr, + m„ (vr, - v^) = 0, 



(4.5^ 



which is the same as the equation of motion of the linear representation (4.49). 



In a similar way, this equivalence can be shown by starting with the linear repre- 
sentation and substituting ni = v{z)ne. 



Pseudoscalar Mass Eigenvalues 

We investigate two ways to calculate the pseudoscalar eigenvalues m^. Rearrang- 
ing and eliminating the longitudinal component is one strategy outlined in [2] 



and is briefly presented here. Alternatively, Appendix A. 3. 2 contains a numerical 
routine we use to solve the set of coupled equations. Using this method, we find 
that TTe, the ratio of tt/ and v{z), is extremely sensitive to boundary conditions. 
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the reason being that v{z) goes to zero as z goes to zero. Resolving the eigen- 
values is difficult and subject to significant numerical error. Fortunately, we have 
shown explicitly that physical results do not depend on the particular represen- 
tation of the pseudoscalars. Therefore we determine the eigenvalues in the linear 
representation. 



First, we follow the method of [2j. We manipulate (4.41) and (4.42) to elimi- 



nate the If field. Adding (441) and (442) yields 



(4.59) 



We then use (4.59) to replace the first term in (|4.42|) and solve for 99, 

1 



-a 



ae 



9^ 2 2 o 



{h{z) ^d^{h{z)dz-Kn)) + mid, Tin 



(4.60) 



where h{z) = a{zYv'^e Using (4.59) again, we denote S^vr vr and rearrange 
to find the eigenvalue equation 



-dz[h ^dz{hiTn)] + gla^v^TTn = mliTn . 



(461) 



This can be put into the Schrodinger-like form by substituting H = n / ^Jh{z). 

-K + V{z)Iln = mlUr,, (4.62) 



Then (4.61) becomes 



where the potential takes the form 

V{Z) : 



3 

4? 



2h z- 

3v' v'^ 30' v'(j)' 0'2 
— + 2— + - + ^ 



zv 



2z 



V 



4 



" + ^.(4.63) 



Z^ 



Solving (4.62) using a standard shooting method gives the mass spectrum shown 



in Table 42 There is no light pseudo-Goldstone boson and no large mass gap 
between the first two eigenvalues. What we have done is taken two second-order 
differential equations and reduced them to a third-order differential equation and 
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found the eigenvalues of S^tt. We do not recover mass values that match the 
experimental mass spectra using this method. 

The numerical routine described in Appendix |A.3.2 calculates the mass eigen- 
values of (4.52) and (4.53), which are then plotted in Figure 4.6 and listed in 



Table |4.2[ Solving the set of equations directly produces a mass spectrum with a 
light pseudo-Nambu-Goldstone boson and mass eigenstates that match well with 
the observed radial pion excitations. These results show that eliminating one of 
the fields from the system of second-order differential equations also eliminates 
information from the mass spectrum. Thus, the eigenvalues of dzTCe do not match 
those of TTi or vTe, at least not in this model. The validity of equating the two sets 
of eigenvalues in [2] may need to be reassessed. 



n 


TT Data (MeV) 


TTi (MeV) 


Large-n tti 


d,7i, (MeV) 


1 


140 


143 




1440 


2 


1300 ± 100 


1557 




1706 


3 


1816 ± 14 


1887 




1925 


4 


2070* 


2095 




2117 


5 


2360* 


2298 


2245 


2290 


6 






2403 


2451 


7 






2551 


2601 



Table 4.2: The observed masses [T] and calculated masses using the linear rep- 
resentations. The large-n limit solutions are valid from n ~ 4. From that point 
onward, the numerical method used becomes increasingly inaccurate and tends to 
skew the ni eigenvalues to larger values than are expected from the linear behav- 
ior. The eigenvalues found using the method of [2] are also shown. *=Appears 
strictly in the further states of [1]. 



For large-ra excitations, the technique described in Appendix A. 3. 2 runs into 



boundary condition problems for n > 4. As the number of oscillations in the 
eigenfunctions increases for higher n modes, the routine finds eigenvalues that are 
skewed to larger values. To uncover the correct asymptotic behavior for large n, 



we take the large- 2; limit of (4.52) and (4.53). As n increases, the eigenfunction is 
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Figure 4.6: The mass spectrum calculated in the AdS/QCD model is plotted 
along with the experimental data pQ. The eigenvalues display two characteristics 
matching the QCD pion spectrum: (1) light ground state and (2) a large gap 
between the ground state and the first excited state. The large-n mass trajectory 
clearly follows our calculated eigenvalues from n ^ 4 when our numerical routine 
inadequately follows the oscillations of the higher eigenfunctions. 



largely determined by the behavior of the effective potential at large z. At large 
z, the dilaton and tachyon behave as 

v{z) = Tz, (4.64) 
(p{z) = fi'^z^. (4.65) 

To take the laige-z limit of both representations, we introduce a new dimensionless 
parameter, ^ = fiz, and expand in C,- In the exponential representation, we find 



that (|4.45|) and (|4.46l) become 



9l 

^2 
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(4.66) 
(4.67) 



Similarly, in the linear representation the expansion of (4.52) and (4.53) at large 
^ yields the eigenvalue equations 



ml 



mlT 



(4.68) 
(4.69) 



where (') indicates differentiation with respect to ^. Each set of equations appears 
to describe a pair of simple harmonic oscillators, the equations of motion of which 



are 



{2k + l)^k. 

{2k + l)TTk fc = 0, 1, 



(4.70) 
(4.71) 



It is a reasonable assumption that yj^ = CfcTr^; it ensures that (4.66), (4.67), (4.68) 



and (4.69) have solutions. Using the form of (4.70) and (4.71) to solve for Ck and 
m\ in both representations, we find 



Ck 



rrii 



9l 



{2k + l)/i2 + gir' 



for the exponential representation and 



Ck 



1^1 



glr^ + {2k + l)f,^' 

{{2k + 3)/i^ - '^kT^) {qIT^ + {2k + l)/i^) 
{2k + l)/i2 



(4.72) 
(4.73) 

(4.74) 
(4.75) 



for the linear representation. 
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So far, we have neglected an important fact: z > 0. The eigenfunctions (p 
and vr describe /ia//-harmonic oscillators and contain only half the modes that full 
harmonic oscillators do; therefore, we must take k — )■ 2k. The mass eigenvalues 
for large n, where n = + 1, in both representations then become 



4,5. 



(4.76) 



which are also listed in Table |4.2 and plotted in Figure 4.6 Combining (4.76) and 
the numerical technique, we obtain all the pseudoscalar eigenvalues. By simple 
investigation, we find that the large-n eigenvalues should be trusted over the ones 
found with the numerical routine for n > 4. 



4.3.3 Vector sector 



The soft-wall model with the quadratic dilaton describes the p meson spectrum 
surprisingly well [18J. In fact, since the scalar field VEV does not couple to the 



vector sector, any dilaton with the behavior (4.4) causes the vector mass spectrum 



to exhibit linear trajectories for the higher resonances. Examining the QCD exper- 
imental data, one sees that the p mass spectrum exhibits linear behavior around 
p(1465) or p(1720); therefore, one expects the appropriately modified dilaton will 
only affect lower lying resonances as higher eigenfunctions localize towards the IR 
and are less dependent on small-2 behavior. 



From the action (4.3), the vector equation of motion is derived exactly as in 
Section 3.1.1 Using the KK decomposition, V^{x,z) = Yl'^=o^"'{z)^^J■<ni^) 



(3.20), we obtain 



1. 



-Vn + miVn = 



(4.77) 



using the axial gauge V5 = 0. Using Appendix A. 2, we transform (4.77) to a 
Schrodinger form 



12 



1 



UJ" I Vr, 



where 



Vn 



UJ 



(z) + log 2:. 



(4.78) 



(4.79) 



Fully simplified, (4.78) becomes 

1 



2z 
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(4.80) 



with the boundary conditions, 



hm Vni^z^) = 
dzV„{z —7- oo) = 0. 



(4.81) 
(4.82) 




Since the dilaton specified in (4.13) is modified for small z, there is a change 
in the slope of the mass spectrum around n = 2 which matches the behavior of 
the experimental data. The numerical vector mass spectrum is compared to the 



experimental data in Figure 4.7 and listed in Table 4.3 While the prediction 
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n 


p experimental (MeV) 


p model (MeV) 


1 


775.5 ± 1 


475 


2 


1282 ±37 


1129 


3 


1465 ± 25 


1429 


4 


1720 ± 20 


1674 


5 


1909 ± 30 


1884 


6 


2149 ± 17 


2072 


7 


2265 ± 40 


2243 



Table 4.3: The experimental and predicted values of the vector meson masses. 



for the p(775) mass is low, the rest of the vector meson masses are in reasonable 
agreement with experiment. Most likely the agreement with the p(775) could be 
improved upon by using a parameterization of v{z) which rises more rapidly to 
its asymptotic value at large z. Nevertheless, the purpose of this analysis is to 
explore the consequences of the improved CSB in soft-wall AdS/QCD models, not 
just to fit data. 



4.3.4 Axial-vector sector 



Unlike the vector field, the axial-vector couples to the scalar field VEV, producing 
a z-dependent mass term in its equation of motion. Similarly to the vector field 
case, the equation of motion assuming A^{x,z) = A'^{x)An{z) using the axial 
gauge ^5 = is given by 



i?2 



-dlAn + uj'd.Ar^ + gl—v\z)An = m\^A. 



(4.83) 



Using the same transformation in Appendix A. 2 as for the vector field, where 
An = e"/^a„, one can express (4.83) as 



-d!a^+[\^"-l^" + 9l^v\z)^ 



V (z) an = rriAan, 



A'2 



1 

2^' 



+ 



+ gl—v\z)]an = ml an. (4.84) 



4^2 
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The expression (4.84) for the axial-vector field matches that of the vector field 
except for the additional term, g1v'^{z)R^ / . Because of this 2;-dependent mass 



term, equation (4.84) is difficult to solve analytically and again requires a numer- 



ical solution. Using the shooting method with the boundary conditions. 



lim an{zo) = 0, 



d.aJz oo) 



0, 



(4.85) 
(4.86) 



the axial-vector meson spectrum is obtained for fixed values of /i, rriq, a, and k. 

The linear behavior of v{z) as z ^ oo leads to a constant shift between the 
vector and axial-vector spectra at high mass values. Comparing the equations of 



motion (4.80) and (4.84) for these fields one finds the asymptotic behavior 



Am^ = ( 



rriA — m 



Vr 



)„- 



d2 

g^—v [z oo) 

Z'^ 



Together with the slope /i^, the mass difference, Am^, indicates that k 



(4.87) 



30, 



although (4.87) is only valid as n — )■ oo. The best global fit to all the data suggests 
K = 15, which is used in this work. The results of our analysis are plotted in Figure 
4.8 and displayed in Table 4.4, The ai(1260) resonance is predicted to within 5% 



and there is good agreement with the higher resonances of ai. 

Note that from (4.87), Am^ > implies that k > 0, which means that the 



potential in (4.3 ) is unbounded from below. To address the stability of the gravity- 



dilaton background requires a complete fiuctuation analysis generalizing the work 
in [37] . Even though this is beyond the scope of the present work it does suggest 
that higher-order terms may be needed for stability. 



4.4 Other Dynamics 

To further confirm the correspondence between our AdS/QCD model and large- 
Nc QCD, we take a closer look at the Gell-Mann-Oakes-Renner (GOR) relation, 
the coupling gp^Trn, and the form factor FT,{q). We find that our model satisfies 
the GOR relation naturally and matches -F^r data quite well, but produces values 
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Figure 4.8: Comparison of the numerical results for the axial-vector mass eigen- 



values using the tanh form (4.14) of v{z) (solid) with the QCD ai mass spectrum 



for gp^TTTT smaller than some soft-wall models and the experimental value. In both 



Sections 4.4.1 and 4.4.2 we need the axial-pion mixing equations that come from 



Section 4.3.2 



dzAp{q,z) A^{q,z) 



glL^v\z) 



0, (4. 



-dzip{z) + 



glR^v\z) 



{ti,{z) - v{z)) = 0, (4.89) 



q^dzViz) + ; ^ dz7i,{z) = 0, (4.90) 
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n 


ai experimental (MeV) 


ai model (MeV) 


1 


1230 ± 40 


1185 


2 


1647 ± 22 


1591 


3 


1930^?^ 


1900 


4 


2096 ± 122 


2101 


5 


2270^^[; 


2279 



Table 4.4: The experimental and predicted values of the axial-vector meson 
masses. 



where we have switched to momentum space and separated out the source from 
the axial-vector field, so that 



AiO,Ro] 

Setting g = 0, we find several facts, 



1. 



A{0,z) = 'f{z)-n{z), 
7[(z) = constant, 



which become useful in our work. 



(4.91) 



(4.92) 
(4.93) 



4.4.1 Gell-Mann— Oakes— Renner Relation 

Using the established equivalence between the exponential and linear representa- 



tions, TTe = 'Ki/v{z), and inserting it into (4.59), we obtain 

-a 



]=mld,<p. (4.94) 
z^ \v / 

Following the method of [15], we construct a perturbative solution in where 
(f){z) = A{0, z) — 1 and use 



ft = ~R- 



(4.95) 



2^0 



From this relation, it follows that 

-k{z) = ml v{z) 



du 



u 
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(4.96) 



The function v?/v'^{u) is significant only at small values of m ~ ^/mja, where 
we may use ( 4.95[ ) to relate the derivative of y4(0,M) to the pion decay constant, 
so that 



2mqCr 



(4.97) 



We find that tt/ = —v{z) solves the axial-vector field equations (4.88), ( 4.89[ ), and 
(4.90) in the region of small z and as g — ?■ 0. As a result, (4.97) becomes the 
expected Gell-Mann-Oakes-Renner (GOR) relation. 



2 £2 



(4.98) 



This perturbative behavior for vTe and tti justifies the use of Neumann and Dirichlet 
boundary conditions, respectively, such that 

7re(0) = -l, ni{0)=v{0) = 0. (4.99) 

The ratio of nig and should be a constant following the GOR relation, 

(4.100) 



J TT 

2^' 



We solve the pair of coupled differential equations (4.52 ) and (4.53 ) for the ground- 
state pseudoscalar mass m^r using differing values of rriq to ensure that the numer- 



ical routine of Appendix A. 3. 2 respects the GOR relation. The results are plotted 
in Figure 4.9 We see a linear relationship between and m^, indicating that 
(4.100) is indeed satisfied. The slope of the line in Figure 4.9 implies ~ 90 
MeV, a result consistent with the input parameters. 



4.4.2 Pion Coupling 



To investigate aspects of the pion coupling, we need the value of /^^ from (4.95). 
For our model where rrig = 9.75 MeV, k = 15, and /i^ = 0.1831 GeV^, we find that 
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4 6 
m^(MeV) 



10 



Figure 4.9: Plot of versus rUg produces a straight line from which the pion 
decay constant is calculated using (4.100). 



/,r = 92.4 MeV. In order to calculate the vector-pion-pion coupling gp„mT, we need 



the Ytttt interaction terms. Using our action (4.3) and the steps and definitions 
in [5], we find that 

where Vn are the KK modes of the p meson. They are normalized as 



(4.101) 



dz 



-Vn{z)Vm{z) = S„ 



(4.102) 



The functions ^{z) and <f{z) must be determined from the system of equations for 



the axial-vector and pion as given in (4.88), (4.89), and (4.90) [5j|. The evaluation 
is made easier because we can approximate by setting ip{z) = ^4^(0, z) — 1 to first 



order. Previous soft-wall models [3] have obtained values a factor of two smaller 
than the experimental result of (^pTrvr ~ 6. Similarly, our calculations produce a 
small value (^ptttt = 2.89. Once we calculate gp^nn, the space-like pion form factor 
can easily be determined from a sum over vector meson poles, 

fn9p 



(4.103) 



where are the decay constants of the vector modes. However, (4.103) converges 
much too slowly and numerically it is much better to use the expression in terms 
of the vector and axial-vector bulk-to-boundary propagators as in [5] 



g^z z^ 



[■k{z) - ^{z))'' 



(4.104) 



The results of our FT^[q ) calculation are plotted in Figure 4.10, and shows a slight 
improvement in matching the experimental values compared to that obtained in 
[5]. It is apparent that the QCD behavior is mimicked reasonably well, beyond 
that expected from a simple soft-wall AdS/QCD model. 



4.5 Conclusion 

In this chapter, we have shown that higher-order terms in the action improves 
upon the soft-wall AdS/QCD model. We correctly incorporate chiral symmetry 
breaking decoupling the sources for spontaneous and explicit breaking. This is 
achieved by introducing a quartic term in the potential for the bulk scalar field 
dual to the quark bilinear operator qq. This changes the dilaton profile for small 
z, while simultaneously maintaining the large z quadratic behavior and therefore 
linear trajectories for the radially excited states. In addition, our model is built 
from the assumption of preserving chiral symmetry for highly excited states, which 
is supported by the experimental values of the QCD mass spectrum. This enables 
us to obtain reasonable agreement within 10% of the QCD meson mass spectra for 
scalar, vector and axial-vector fields, although the lowest lying p and /o predictions 
are not as good. 
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Figure 4.10: The line shows the predicted space-hke behavior of the pion form 
factor F-,^{q^) which is compared to the experimental data obtained from [5j. The 
triangles are data from DESY, reanalyzed by |6j. The diamonds are data from 
Jefferson Lab The circles [H] as well as the star [H] are also data obtained from 
DESY. 



Even though our modification of the soft- wall version of the AdS/QCD model 
is simple and predictive, any further progress must recognize the limitations of this 
type of phenomeno logical model. Genuine stringy behavior is most likely required 
to fully describe the characteristics of QCD. Nevertheless, some features such as 
masses and couplings seem to agree better than expected and it would be worth 
using the modified dilaton profile to study further details of the meson spectrum. 
On the theoretical side, it would be interesting to further understand the soft-wall 
model from the top-down including finding a dynamical solution of the features 
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exhibited in our model along the lines considered in [HS]- In addition, the stability 
of the scalar potential will most likely require higher-order terms that can only 
be studied from the top-down. It is interesting that the simple five-dimensional 
model contains QCD-like features and suggests that a further understanding of 
QCD can be obtained from the gauge/gravity correspondence. 



Chapter 5 



Finite Temperature 
Thermodynamics 

"Nothing in life is certain except death, taxes and the second law of 

thermodynamics . " 
-Seth Lloyd 

5.1 Introduction 

In this chapter, we focus on two descriptions of our QCD-like gauge theory. Con- 
fined hadronic matter where baryons and mesons are the primary degrees of free- 
dom is characterized by the AdS metric similar to the one from Chapter |4j De- 
confined plasma where quarks and gluons associate freely is modeled by extending 
the metric to an AdS-Schwarzschild solution, describing an extra-dimensional gen- 
eralization of a black hole. The thermodynamic description of black holes is well 
established [751 EHl EZ] • The correspondence connects the black-hole thermody- 
namics with the thermodynamic properties of a deconfined gauge theory. The 
transition between these two phases occurs at a critical temperature Tc and corre- 
sponds to the Hawking-Page transition in the gravity dual. Experimental evidence 
confirms that the QGP is strongly coupled [IHl [13, UHl [IS] , lending support to the 
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theoretical picture obtained via a five-dimensional dual gravity description. 

In this chapter, we find a finite-temperature description of a strongly coupled 
gauge theory with a soft- wall geometry based on the model in ^5]. The five- 
dimensional Einstein equations involve two scalar fields. The black-hole metric 
is asymptotically AdS with an event horizon located at a finite value of the bulk 
coordinate, z. In the string frame, one scalar field plays the role of a string- 
theory dilaton, while the other scalar field is similar to a string-theory tachyon. 
In terms of the correspondence, we begin with the assumption that the dilaton 
field is dual to the dimension-4 gluonic operator, "YxiF^yF^^^ while the second 
tachyon field is dual to the dimension-3 chiral operator qq. The true field/operator 
correspondence is more nuanced, as we will discuss. 

The solution outlined in this chapter is an expansion only valid in the limit 
of Zh < 1. Unfortunately, this is the price we pay to study the soft-wall ther- 
modynamics; the equations of motion have no known closed-form solution nor 
a reasonable numerical solution. We explore a purely phenomenological model 
that does not dynamically generate the soft-wall geometry in Appendix |C.3 We 
attempt to introduce rigor into the soft-wall thermodynamics, beyond what has 
been done in [8ntlTn2]. 

The thermodynamics resulting from our solution leads to interesting conse- 
quences. We show that a nonzero thermal condensate function induces a 
phase transition [HH ESI E] • Furthermore, Q contributes leading-order terms that 
are absent in the lattice results to the soft-wall thermodynamics |103l 1104] . How- 
ever, qualitatively the two agree. The entropy has the expected behavior at high 
temperatures, scaling as T^, and the speed of sound through the thermal plasma 
is consistent with the conformal value of 1/3, agreeing with the upper bound 
advocated in |105j . 



Our analysis begins in Section |5.2[ where we present the thermal AdS and 

we 



black-hole AdS solution and compute the on-shell action. In Section 5.3 



study the thermodynamics of our solution, including a general expression for the 
entropy and squared speed of sound. We then calculate the free energy difference 
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by carefully matching our two solutions at the AdS boundary. This enables us to 
compute the transition temperature between the confined and deconfined phases 



of the gauge theory. Concluding remarks are given in Section 5.5 



5.2 Finite- Temperature Action 

We begin by specifying the five-dimensional action in the string frame and trans- 
forming it to the Einstein frame. The equations of motion become simpler in the 
Einstein frame. We are extending the AdS solution with two scalar fields first 
investigated by |95j. 



5.2.1 5D Lagrangian 

We start with the string-frame action inspired by the dimensionally-reduced type 



IIB supergravity action first introduced in Section 3.2.1 



1 



string 



167rGR 



d Xy/—g 



e - I R + Ag^'^dM^dN^--g^'^dMxdNX 



(5.1) 



where £meson contains all the mass terms not considered in this chapter, is 
the string-frame scalar potential, and the indices M, = (t, xi, X2, 0:3, z). We are 
working with two scalar fields $ and x where previous models have only focused on 
the dilaton |671 [H |85l EHl |89l EB [IO6] . The inclusion of a Gibbons-Hawking term 
Sq\j is an attempt to be more rigorous than previous models. The string-frame 
metric is assumed to have an AdS-Schwarzschild form 



ds = QMNdx dx 



^ ^ - j{z)dt^ + dx^ + '^^ 



(5.2) 



where f{z^ determines the location of the black-hole horizon. Furthermore, the 



string-frame action (5.1) contains similarities with noncritical string theory. The 
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$ scalar field behaves like a dilaton, while the x scalar field behaves like a closed- 
string tachyon field. This suggests that our setup could be the low-energy limit 
of some underlying string theory. 

While the string- frame action provides a suitable starting point, it is more 
practical to do calculations in the Einstein frame. Switching to the Einstein 
frame involves a simple conformal transformation, 

9mn ^ ^^^9\A}- (5.3) 
The gravity-dilaton-tachyon action in the Einstein frame then becomes 

Se = j d'x^g [r - \{dct>f - \{dxf - V{ct>, x)) + Sgh. (5.4) 

where (f) — •\/8/3$ and V — Vget*. We explicitly define the Gibbons-Hawking 
term Sqh as 



d'x / dT^/^K, (5.5) 



457r2i?3 

where K — j'"^K^iy and 7 is the four-dimensional induced metric at the AdS 
boundary. The extrinsic curvature K^i, is defined by 

= V^n^ = ^n'^dulnu , (5.6) 
where the vector hm is the outward directed normal to the boundary, normalized 

by 

^MivnV = l. (5.7) 

The boundary of the AdSs space considered here is the z — plane, meaning the 
normal vector is 

^/g7z \ozJ 



The Gibbons-Hawking term does not affect the equations of motion, but will have 
consequences when considering the free energy and deconfinement temperature. 



To convert to a description at finite temperature, we compactify the Euclidean 
time coordinate, t = it ^ it + where /3 = 1/T is the inverse temperature. The 
finite-temperature metric in the Einstein frame then becomes 



a{zy f{z)dT' + dr + 



(5.9) 



with the finite-temperature action given by 
1 



IGttG, 



d X 



1 



8nG. 



dr 
d^x 



dzy 



■{dxf-V{ct),x) 



dTy/—'-)K. 



(5.10) 







The Einstein frame metric (5.9) and action (5.10) will be used to solve the 
Einstein's equations in two realms, thermal AdS (thAdS) and black-hole AdS 
(bhAdS). All quantities with a subscript o are associated with the thAdS solu- 
tion. 



5.2.2 Thermal AdS Solution 

We begin by obtaining a thAdS solution that is effectively equivalent to the zero- 
temperature case explored in Chapter |4| where f{z) = 1. Assuming the scalar 
fields are a function of only the z coordinate, the action can be expressed as 

+ S^,GH. (5.11) 

where V is the spatial three-volume and a UV cutoff at z = 5 has been introduced 
to regularize any singular behavior at the AdS boundary. 
The metric associated with the thAdS solution is 

ds" = a{zf{dT^ + dx"" + dz^) = e-2^'^(^)^(rfr2 + dx^ + dz^), (5.12) 

z^ 
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where the constant c depends on the conformal transformation. In our case, 
c = 1/a/6. Two equations come from the five-dimensional Einsteins equations, 

,/2 



a a 
12— -Q — 
a 

a'2 a" 
6^ + 3 — 



-a^V{z) 



a 

and two more equations come from the varying the action. 



,dV 



a 



,dV 



x" + 3x' 



(5.13) 
(5.14) 

(5.15) 
(5.16) 



where prime (') denotes derivatives with respect to z. Using (5.12), we express 



(5.13), (5.14), (5.15), and (5.16) all in terms of 0, Xi and V {(f){z), xi^))} 



X'' 
V{z) 
dV 

dV 
dx 



z 



+ 



2 A" 
Z (p 



z^e^'^ 



i?2 



X 



T z 



(5.17) 
(5.18) 

(5.19) 

(5.20) 



where we see that the nonlinear term 0'^ has conveniently cancelled in (5.17). Of 
the four equations, we find that three are independent. 

As shown in [95], the solution in the soft-wall model with a quadratic dilaton 



gives 



2 2 
[I Z , 



X{z) 
V{z) 



3 

2^[iz 



(-12-20/i2z2_i6^V) 



(5.21) 
(5.22) 
(5.23) 
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where n sets the hadronic mass scale. The quadratic behavior of the (f) solution 
(5.21) leads to a Regge-like hadron mass spectrum (m^ ~ /i^n). The linear be- 



havior of X indicates similar chiral symmetry consequences as the scalar VEV 
function v{z) did in Chapter |4j The potential is a function of the z coordinate 
with no unique solution for V{(j),x)- However, a potential was suggested in 



|eT^ + 202e^'^ - 12 (^3e^ - 2 |^1 - -^^ e^^ 



(5.24) 



We can define an alternative potential that still produces (5.21), (5.22), and (5.23) 



2 

e^ 



7 



(5.25) 



3v^ 27' 

More examples of consistent and well-defined potentials written in terms of two 
scalar fields can be found in |107] . 

Part of the the free energy expression is found by substituting the thAdS 



solution into (|5.11|) and finding the Gibbons-Hawking term. In general, we find 

8a" 4a'2 



the Ricci scalar is 



R 



(5.26) 



For the thermal AdS solution, we find the extrinsic curvature takes the form of 



1 6 



M 



Calculating the trace, we find 



(5.27) 



(5.28) 



Substituting this value into (5.5) gives the on-shell Gibbons- Hawking term for the 
thermal AdS solution, 

V 



So. 



GH 



457r2 R^T 



(8ao«o) • 



The on-shell action (5.11) then becomes 

V 



457r2 R^T 



GH 



NlV2a'{5)a'{5) 



157r2 



R^T 



(5.29) 



(5.30) 
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The on-shell action is a pure boundary term and strictly depends on the AdS 
boundary conditions, where 6^0. 



Field/Operator Correspondence 

According to the AdS / CFT dictionary, the scalar field masses determine the cor- 



responding operator dimensions in the dual gauge theory via (2.18). Expanding 
the potential (5.24) from [93], we see that 



suggesting that (p is dual to a dimension-2 operator, and x is dual to a dimension- 



3 operator. When compared to the standard form (2.16), the behavior of the 



solutions (5.21) and (5.22) confirm the operator dimensions suggested by (5.31). 



However, the potential (5.25) gives the masses as 

" " 0, 



(5.32) 



indicating that (p and x ^ire dual to a dimension-4 and dimension-3 operator, 
respectively. It is fairly clear that the chiral operator, qq, and the gluonic opera- 
tor, Tr[F^], are the dimension-3 and dimension-4 operators, but the dimension-2 
operator is much less clear. 

The most likely dimension-2 operator candidate is A^, which becomes a lo- 
cal expression in the Laudau gauge, d^A^ = 0. Coupling a source term to Aj^, 
however, makes the theory nonrenormalizable at the quantum level. A quadratic 
source term can be added to remedy this obstacle, though, this ruins the energy 
interpretation of the effective action [108]. In the context of the AdS/CFT cor- 
respondence, Aj^ is often understood to convey information about the topological 
defects in the gravity dual |109l 1110] . Much more work concerning Aj^ has been 
conducted in [ml fTT2|[TT3]. 

In the current soft- wall case, the field/operator correspondence appears compli- 
cated. The ambiguity stems from the fact that the original AdS/CFT dictionary 
was formulated considering purely free scalar fields. The potentials (5.24) and 
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(5.25) clearly have interaction terms. Resolving the issue of whether interaction 
terms affect the field/operator correspondence and determining the interpretation 
of the dimension-2 operator is ultimately beyond the scope of this thesis. In that 



spirit, we expand upon the published potential (5.24) and assume that the fields 
X and correspond to the operators qq and Tr[F^], which is the general consen- 
sus. The significance of our work relies on the fact that the dilaton is dual to 
some temperature-dependent operator. The identity of that operator is a topic 
for further research. 



5.2.3 Black-Hole AdS solution 



Next, we consider the black-hole AdS solution that describes a deconfined phase 
mimicking a free gluon plasma. Assuming the solutions are only a function of 
the z coordinate, the five-dimensional action associated with the black-hole AdS 
solution simplifies to 



V 



457r2 R^T{zh) 



dz 

+ Sbh,GH, 



(5.33) 



where Zh is the location of the black-hole horizon. We will see that Zh is directly 
related to the temperature of the gauge theory. We begin with the black-hole 



metric (5.9) and find four independent equations, 

a'{z) 



,dV 



-3/'(^) 



a{z) 



12 



a{zy 
a{zy 



d(j) 
,dV 

dx 



a\zf a"{z) 

o 

a{zY a{z) 

f{z)f{z) + f{zy 



[z) + 3f{z 
f{z)x"{z) + f{z)x'{z) + 3f{z)x 



a{z) ' 
,a\z)_ 
a{z) 



(5.34) 
(5.35) 
(5.36) 
(5.37) 



Unlike in the thAdS case, the potential is already determined. We must use the 



potential (5.24) to connect the soft-wall action to the free energy investigated in 
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Section 5.3 With four independent equations and four unknown functions, f{z) 



a{z), (f){z), and x{^)y "we no longer can assume a fixed relation between the warp 
factor a{z) and the dilaton 4>{z). These quantities must evolve independently 
as the temperature varies. The system of equations associated with bhAdS are 
difficult to solve but for the simplest cases. The conformal case discussed in 
Chapter |3] is the only known exact solution. 

Using the series expansions, we construct another solution. We use the thAdS 



solution of Section 5.2.2 as the starting points for these series expansions, 



aiz) 



z ' 

oo 



n=2 

oo 

n=l 

oo 

i + J2fn{g,z,)z^, 



(5.38) 
(5.39) 

(5.40) 

(5.41) 



n=4 



where we see that one of the black-hole conditions, /(O) = 1, is automatically 
satisfied. The condensate function Q{zh) plays an important role in the free energy 



and phase transition of the thermal plasma. By solving (5.34), (5.35), (5.36), and 



(5.37) in successive powers of z, we find the coefficients up to n = 8 in terms of 
fiiQ, Zh)- the other black-hole condition, f{zh) = 0, determines the final unknown 
coefficient. We construct a solution where the majority of the coefficients are zero. 
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We calculate the non-zero coefficients for the metric, 



y/6 



m2 = -iT^Gizh), (5.42) 



^(^OM, (5.43) 



= 8;.^ + 3V6^(..) I, - + 2l7I^^^^'^^ ^ lisTT^^^'^^ 

71 x2^ 2543/ , 1507 , 11773 ^ 3 

+^/4^(..) + ^^(-.) + + ^^^JM^^) 

^^^'^'^^'-^s^^^^ri (5.44) 
the field (p, 

Pa = -G(zh), (5.45) 

P6 = -4/4-^^(^/.) + ;r4i^e^(^/^)'> (5-46) 



1 / 80vV, 12/x\^. . 1408/iV , , 

+^fMzH) - -^S{z,) - -g{z,) 

M-^f\^ (5-47) 



100876/^4 
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the field x, 



and f{z) 



ci = ^/2V + ^g;(^,), (5.48) 

C3 = . ^ (5.49) 

^8/.2 + ^^(^,) 

V + 3v^6^(^/.) V 

8/i4 + 3y6^(z,) V 14^2 189 

lllv^.^. n2 1123/i 2 67 3 

15iyj 
22V 



+ (5.51) 



/4 — — 7 7 \ (5.52) 

h = ^U+-^J,g{z,), (5.53) 



/8 = fJ^'h + ^fM^k) + ^JMzhY. (5.54) 



In general, an arbitrary power z"- multiplying a factor of g{zh) in the coeffi- 
cients Pn and rrin has direct consequences for the free energy. For < n < 4, 
the free energy contains divergences of the power n, while for n > 4 the function 
g{zh) does not affect the free energy. Only the n = 4 coefficients affect the free 



energy expression. As we will see in Section 5.3, the condensate function Q then 
plays a crucial role in giving rise to a finite transition temperature. 
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The calculation of the Gibbons-Hawking term is more complicated in the 
bhAdS solution. The induced metric in this case is 



^ = a{zf {f{z)dT^ + dx^) . 



(5.55) 



To obtain the trace of the extrinsic curvature, we need to separate the r compo- 
nents from the rest, 



Using (5.27), we find that 



00 



2,/Ja 



(5.56) 

(5.57) 
(5.58) 



Given that = a^\fj ^ the Gibbons-Hawking term becomes 



S> 



bh,GH 



(8/ a^a' + / a^) . 



(5.59) 



Finally, an expression for the on-shell action can be obtained by substituting 



the solutions back into ( 5.33 ) . Using the fact that the Ricci tensor can be expressed 

as 

4/a'2 8fa" 8f'a' f" 
R=-—. V-^, (5.60) 



the on-shell action can then be written as 

N^. V 



V 



,GH 



457r2 R^T 



{6fiS)a'{S)a\S) + f'iS)a%6)) . 



(5.61) 



Again, the on-shell action is a pure boundary term and depends only on the AdS 
boundary conditions. This action will be used to compute the free energy in 
Section 15.31 
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Field/Operator Correspondence at Finite Temperature 



The scalar field solutions (5.39) and (5.40) correspond to turning on thermal 
condensates in the gauge theory. To show the relation between Q{zh) and the 
operator condensates, we assume the kinetic term of a canonically normalized 
bulk scalar field fluctuation, denoted u, to be 

5bulk = ^ / d'xy/^ dMCod^cu. (5.62) 

Let us first consider the scalar field (p, where the coupling to the four-dimensional 
boundary operator is 

^boundary = y" C^^a; W0 TrF^ (5.63) 

In terms of the 't Hooft coupling A = e*, the Yang-Mills field strength is 

^boundary = - J d'^X ^TrF^, (5.64) 

making the dilaton fluctuation, 

55boundary = \j d^x 5$ e-*Tr(F2). (5.65) 

We are only interested in the fluctuation 5$ = $— $0; which allows one to compute 
the difference between thermal and vacuum values of (TrF^) [84]. From (5.39), we 
find that 5$ = —^J'i/9>Q{zh)z'^. Recall that the relation between the expectation 
value of a A-dimensional operator in d-dimensional space and the field is 

4-4^- (5.66) 

The function Q{zh) then relates directly to the gluon condensate, 

{T.{F^f) - {Tt{F^)\ = -yl (^H^) Oiz,). (5.67) 

A similar correspondence can be obtained for the scalar field Xy where it is 
tempting to relate x to the three-dimensional operator qq. In this case, the coef- 



ficient of the z term in (5.40) corresponds to the renormalized chiral condensate. 
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{qq) — (gg)o- For this correspondence to hold, x would need to be a bifundamental 
field in the gauge theory. This can be done by promoting x to a bifundamental 
field x"*; where a, h are group indices and writing, for example 

{X''') = x{z)\ : : I (5.68) 




for an SU{2) x SU{2) symmetry. The bulk x solution (5.40) is then identified 



with the vacuum expectation value xi^) in (5.68). 

Under this assumption, we can follow the same procedure as for the field, 
and derive an operator correspondence for x- Assuming a bulk kinetic term of the 



form (5.62) and boundary operator coupling 

<Sboundary = J d'^X U^qq, (5.69) 

we find that 

m-m.^-I^C-^-l'^). (5.70) 



277r2 \ 12 8 /i5 ^ 
Thus, we see that the function Q{zh) is intimately related to the thermal conden- 
sates (qq) and (TrF^). In fact, the ratio of the gluon condensate and the chiral 
condensate is given roughly by 

(Tr(F")2) - (Tr(F")2)o ^ 16 

— — ^ — V6A/i, 5.71 

{qq} - {qq)Q 5 

and is consistent with that obtained in perturbation theory to first order |114j . 
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5.3 Thermodynamics 

Having determined the black-hole solution of the five-dimensional gravity-dilaton- 
tachyon system we now investigate the thermodynamics. We begin with Hawk- 
ing's black-hole thermodynamics. An expression for the temperature T was de- 
rived in Section 13.21 

djjz) 1 



T(zh) 




2 "P ^/i 4v^m2 12fj.2+8fiizl+3V6zlgizh) . 

Of course, all thermodynamic quantities depend on the condensate function, 
whose behavior we address in Section [5^ 

As stated earlier, the graph of T{zh) contains information about the stability of 
the black-hole solution. The specific heat ch of the system is found using standard 
thermodynamics, 

dE dE/dzh , . 

'^=dT = djW. ^ 
In general, E'{zh) < 0; therefore, the zeros of T'{zh) indicate where the specific 
heat switches sign. Whenever T'{zh) < 0, we find a stable AdS black hole. The 
behavior of the temperature as a function of z^ with and without a condensate 



function is shown in Figure 5.2 We use (3.151 ) to indicate the general the large- 2;/, 
behavior. 



5.3.1 Entropy 

The entropy is found in the usual way from black-hole thermodynamics, a subject 
extensively covered in [7S1 [THl E3- expected for a relativistic gas at high 
temperature, the entropy density, s(T), behaves like T^. In our model, we can 
also calculate the subleading temperature behavior and check that it is consistent. 
We compute the entropy density from the area of the black-hole horizon using the 
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induced metric 7. Using the metric (5.38), we obtain 



A 



bh 



where V is the spatial volume. Using the first few terms in our metric expansion, 
the entropy density is then given by 



5.3.2 Speed of Sound 



With the temperature and entropy known, we also determine the speed of sound 
through the deconfined medium of the gauge theory. The speed of sound Vs char- 
acterizes the hydrodynamic evolution of the deconfined, strongly coupled plasma. 
It has been suggested that vl in holographic models obeys an upper limit of 1/3 
[105] . This can be checked for our solution using the relation 



V 



2 ^ ^ d\ogT 

dlogs 



3 a log a 

In our model, the exact form of v"^ greatly depends on the form of Q. However, 
we do confirm that 

- -4, (5.77) 

in the high temperature limit, recovering the conformal limit with an upper bound 
of = 1/3. 

5.3.3 Free Energy 

The free energy of the deconfined phase is calculated in two ways. We can use 
thermodynamic identities to define the free energy as 

T=- J SdT = T^in- J s V^^dzh. (5.78) 
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where we have found the Bekenstein entropy and have an expression for T{zh). We 
have exphcitly written the integration constant J-min since the problem with using 



(5.78) is common among energy definitions, setting the zero point. Because Zh{T) 
is a multi-valued function, the integral is non-trivial. Calculating the value J^uim 
generally requires the expression of entropy in the large- z/j region, exactly where 



our expanded solutions are invalid. While the integral in (5.78) is completely 



calculable, we find the free energy by using the on-shell action. Finding the zero- 
point energy using actions is much easier; one merely subtracts the background 
free energy as defined by the thAdS action, 



J^ = \ims^oT[S,h{S)-So{S)]. 
Both Sbh and Sq have been computed earlier. 



(5.79) 



Before evaluating (5.79), we must properly match the thermal AdS and black- 



hole AdS metrics at the boundary, 5 — ?■ [93l |8l]. This requires matching the 
intrinsic geometry of the two solutions at the boundary cut-off, where 



ao(5) = a(5)v^, 
Voao{6f = Va{6f. 



(5.80) 



In order for (5.80 ) to be satisfied, we must evaluate the thAdS and bhAdS solutions 



at different cut-off points, 6 and 6 respectively. We find that 



52 



(5.81) 



Combining the matching with (5.79), we obtain a rather simple expression for the 
free energy. 



J" 



457r2i?3 



lim 



(5^0 



6f{S)a{Sra'{S) + f'a{Sr - eVWHS)'^ 

ao(5)2 



which can be reduced to 



457r2 



r- \ 2N'^V 



-TS. 



(5.82) 



(5.83) 
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It should be noted that we also checked that the black-hole energy E satisfies the 
thermodynamic formula E = T + TS by computing the ADM energy in Appendix 



5.4 The Condensate Function 



All the thermodynamic relationships rely on the behavior of the condensate func- 
tion; therefore, we need to find a solution to Q to evaluate the temperature, 
entropy, and speed of sound. The free energy of the system gives us enough infor- 



mation to solve for Q. We only need to set (5.78) and (5.83) equal to one another. 



Taking the derivative with respect to removes any unknown constants, giving 

dT 



dzh dzh 



sV 



dzh 



or in more simplified terms, 
dQ{zh) 



where 



dzh 2V6V{ 



V{z) 



Zh) 



a{zh) 



dx a{xy 



^V{z^)) 



(5.84) 



(5.85) 



(5.86) 



Unfortunately, (5.85) is a stiff equation. Stiff differential equations contain terms 



that lead to rapid variations in its solution, causing numerical instabilities. One 
can often find stable solutions to these equations within a certain region, but our 
case is barely within a region of stability. 

We use two methods to solve for Q. First, we introduce a series expansion for 



G{zh) = gjfi- 



hi 



(5.87) 
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and expand (5.85) in powers of Zh- Performing the expansion and matching coef- 
ficients, we find that the nonzero g/s are 



g-2 = g = 1.43290, 
go = -6.09417, 
c/2 = -29.7867, 
= -570.637. 



(5. 



Since the nonhnear nature of (5.85) spoils the series solution quite quickly, we 



find that only the first two terms give an accurate solution for G in the range of 
Zh < 1, 



G{zh) = 1.43290^ - 6.09417/. 



(5.89) 



Using the leading order terms to inform the boundary conditions, we then numer- 



ically solve for Q in (5.85). Numerically, (5.85) is difficult to solve because of its 



stiff nature. However, we are able to obtain a valid solution in a small range of 
Zh, which matches the series expansion, until Q diverges. The two solutions are 



plotted in Figure |5.1[ With a solution to Q, we can take a second look at the 
thermodynamics. 



Given the expression (5.41), the temperature can be written as 
2 /'l6 + 12^6^ + 27^2 



nzh) 



P Zh 



TTZh \ 32 + 12^6^ + 27^2 
'2048 



7r(32 + 12^6^ + 27^2)2 
+ 516^6^ + 192/ + 256^6^0 + 1728^^0 



+ 216^6^^ 



^0 



(5.90) 



We clearly see that the condensate function contributes a finite piece to the tem- 
perature expression. Increasing the strength of the condensate produces 

2 



lim T{zh) 

g-5-oo nZh 



(5.91) 
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Figure 5.1: The series and numerical solution for Q is plotted. We use the series 
solution for the thermodynamics since the differential equation clearly gives a 
numerically unstable solution for z > 0.3 GeV~^. 



which is a factor of 2 larger than the conformal limit in (3.150) or that found in 



We find that other work in this area has assumed that the condensate 
terms are suppressed logarithmically. However, in the construction of this model, 
the soft-wall set-up has no natural means of generating the logarithmic suppression 
in ^. As a result, we see that the condensate contributes to leading order behavior 
in the small- z/i/large-T limit. 

It will be useful to have an inverted function Zh{T) which we use to transform 
functions of Zh to functions of T. For simplicity, we use only the first-order term, 

irT I 32+ 12v^9 + 2792 j 
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Figure 5.2: The temperature as a function of the horizon location Zh with and 
without a condensate Q. The inclusion of a condensate in this particular model 
precludes a second, unstable black-hole solution from developing. 



Analytically, the entropy of the deconfined phase can be written in terms of 



Zh, 



s{zh) 



1 



457r (32 + IQy/Qg + 27c/2)2 
3623 , 



86^6 279 2 313\/3 

32 q H q ^ 

7 ^ 14^ 14v^ 



+ 



168 " 224v^ 
1.77 25.35/i2 



313 , 90081 
3584 ■ 



9-^ + 



+ 



Zh 



(5.93) 
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or in terms of T, 



AN'^V TT^T^ ( , 17664^6 , 43200 ^ 



s{T) ^ —2 4096+ ;r^9 + 

457r (16 + 12^6^ + 27c/2)3 V 7 

3056\/6 o 57520 . 22331^3 , 102549 ^ 



7 ^ 21 ^ 14v^ ^ 16 

992899^3 . o 9652689^3 „ 65669049 m', 
+ ^9 + 12069/ + + + 



112v^ 1792^2 28672 

« 13.88T^- 50.35/i2T + ..., 

s(T) = C(T)T3, (5.94) 

where the function C(T) modifies the entropy behavior at low temperatures. Our 
result, plotted in Figure |5.3[ agrees qualitatively with the lattice data presented 
in |103] . The high-temperature limit of s(T)/T^, however, is shifted from lattice 
results because of the contribution from the condensate function. 

The speed of sound through the QGP-like thermal phase can be expressed as 

vl{zh) ^ ^ - 2.718/^24 + 21. 938/iSt (5-95) 

vl{T) ^ 1 - 0.689^ + 1.409^. (5.96) 

We clearly see that (5.95) and (5.96) give the expected = 1/3 in the small-^ 



and large-T limit. The speed of sound through our QGP is plotted in Figure [5^ 
As the figure shows, f f deceases more rapidly with the condensate terms included; 
however, the difference between the two cases are not as stark as it was in the 
entropy case. Again, our results would match those of the no condensate scenario 
if the terms contributed by the condensate were logarithmically suppressed as in 

IMI. 



An issue arises when we consider the behavior of Q. In Section 5.2.2 , we argued 
that the dual operator corresponding to the dilaton field is either a dimension-2 or 
a dimension-4 gluonic operator. On dimensional grounds, previous work asserts 
that ^ ~ T^. However, our model clearly shows that any self-consistent solution 
in the soft-wall model requires that ^ ~ T^. The consequences of the temperature 
dependence on the field/operator duality are subjects for further research. 
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Figure 5.3: The entropy density divided by as a function of temperature. At 
high temperatures the entropy density slowly evolves to the conformal case of 
s ~ in each case. 



5.4.1 Deconfinement Temperature 

Finding the deconfinement temperature requires examining the free energy of the 



deconfined phase in (5.83). The point at which J-" = occurs when the black- 



hole solution becomes more energetically favored than the thermal solution. It is 
instructive to first consider the case without a condensate. By setting Q{zh) = 



in (5.83), the transition temperature is determined from the condition 



fi{Q{zh),Zh) = 0. 



(5.97) 



However, (5.97) has only one possible solution, — )■ oo; therefore, no transition 
temperature exists since z/j — )■ oo occurs in the region of the unstable black hole. 
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Figure 5.4: The squared speed of sound in the strongly coupled plasma as a 
function of the temperature. The conformal limit of f ^ = | is recovered at high 
temperatures. 



This mimics the scenario considered in [Sni 1102] . but we include the full back- 
reaction of the scalar field and take the Gibbons- Hawking term into account. 
Therefore, we see that a nonzero condensate is needed to obtain a transition 
temperature in our model. 

Using the leading behavior and the numerical solution, we find the free energy 



behavior and plot it in Figure |5.5[ The phase transition occurs at a critical 
Zc = 0.5262. Using (5.72) and the first two terms of the expansion for Q, this 
corresponds to a critical temperature of = 919 MeV. This is much larger than 



either theoretical reasoning or lattice calculations suggest |114l I115[ I116[ I103[ 1104] . 
There are no current plans at the LHC for experiments that reach such high 
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temperatures; however, experimental evidence for QGP has aheady been detected. 
The large transition temperature is most likely the product of the crude soft- 
wall model and the absence of any suppression terms in Q. As Zh increases (T 
decreases), we see that free energy goes to zero as in [81]; however, we should only 
trust the free energy up to the point of validity of the bhAdS solutions, ^ 1. 
Thus, the true behavior of the free energy at increasing values of Zh requires more 
numerical work. 




Figure 5.5: The free energy J-" plotted as a function of Zh- Only when including a 
gluon condensate do we find a solution to J-" = 0. The free energy with no gluon 
condensate approaches zero as z^ ^ oo, but never crosses the x-axis. 
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5.5 Conclusion 

We have presented a five-dimensional, gravity-dilaton-tachyon, black-hole, dynam- 
ical solution that represents a dual description of a strongly coupled gauge theory 
at finite temperature. The solution generalizes the soft-wall geometry considered 
in [95], which generates a quadratic dilaton and Regge mass trajectory (m„ ~ V^). 
In the string frame, the dilaton and tachyon fields appear to be duals to a gluon 
and chiral operators, respectively. However, the actual field/operator correspon- 
dence is ambiguous. The black-hole solution describes a deconfined free-gluon 
phase. A transition from thAdS exists, predicted at a extremely high Tc = 919 
MeV. A condensate function is needed for any transition to occur. 

Despite being a mere series expansion, we argue that it gives insight into 
thermodynamics of a strongly coupled gauge theory produced by the soft-wall 
model. We expected to find two major contributions to the thermodynamics: 
the underlying conformal limit and condensate terms. Although we considered a 
finite number of terms in our infinite solution expansion, our calculations suggest 



that the condensate terms produce convergent quantities, as shown in (5.91). 
With further study, we believe a viable closed-form, black-hole solution with a 
lower transition temperature will be found. A much more reasonable transition 



temperature is obtained by the phenomenological model in Appendix C.3 

While our model has some features reminiscent of QCD at finite temperature, 
it still represents a crude approximation with a number of shortcomings. The 
bhAdS solution is only valid in the region of small Zh- The power-law dependence 
of the scalar fields and metric does not include the logarithmic corrections needed 
to suppress condensate terms as in |S1], resulting in noticeable changes in the 
behaviors of the temperature, entropy, and speed of sound. Furthermore, using 
the AdS/CFT dictionary, the dilaton field appears to be dual to a dimension-2 
operator, in contradiction to the standard assumption of a dimension- 

4 operator, Tr[(F")^]. Nevertheless, our five-dimensional dynamical black-hole 
solution with two scalar fields provides a toy model to understand the nontrivial 



properties of strong-coupled gauge theories at finite temperature. 



Chapter 6 



Discussion and Conclusion 

"I've learned that people will forget what you said, people will forget what you 
did, but people will never forget how you made them feel." 

-Maya Angelou 

We began this thesis with a quick overview of gauge / gravity dualities. In Chap- 
ter [2| we reviewed basic concepts of AdS/CFT and its connection to AdS/QCD. 
We took a new perspective on the mass spectra produced by the hard-wall and 
soft-wall AdS/QCD correspondence in Chapter |3} In the same chapter, we in- 
troduced thermodynamics into the correspondence via a black-hole metric. In 
Chapter |4| we explored the consequences of adding higher-order interaction terms 
in the action. We found a new way of expressing spontaneous and explicit CSB. 
The modified background dilaton field produces meson mass spectra that agree 
with experimental results better than the previous pure linear trajectories. We 
also edified the pseudoscalar sector that had been neglected by past soft-wall work. 
In Chapter [5} we investigated the means of dynamically generating the soft-wall 
from the gravitational action. We explore an approximated series solution that 
produces a quadratic dilaton and linear tachyon to leading order. 

The phenomenological soft-wall AdS/QCD correspondence mimics QCD in 
a number of aspects: (i) radial Regge trajectories, (ii) running coupling con- 
stant, (iii) Gell-Mann-Renner-Oakes relation, (iv) deconfinement temperature. 
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and (v) high-temperature thermodynamics. In our research, we have shown that 
non-trivial improvements are made when adding higher-order terms. Low-mass 
resonances that experimentally show non-Regge behavior are reproduced. Inde- 
pendent sources of spontaneous and explicit chiral symmetry naturally arise. The 
thermodynamics of the soft-wall model is much more nuanced. Solving Einstein's 
equations that dynamically generate the soft-wall has proven exceedingly diffi- 
cult. Predictions from such attempts give the right high-temperature behaviors, 
and a transition between confined matter and deconfined plasma occurs. How- 
ever, results are skewed by leading order contributions from condensates that are 
thought to be logarithmically suppressed. Phenomenological soft-wall models ap- 
pear to describe the thermodynamics of deconfined QCD much better, as shown 



in Appendix C.3 



Perennial limitations and shortcomings plague some aspects of the soft-wall 
AdS/QCD correspondence. Currently, the running coupling of QCD cannot be 
analytically calculated at low energy. Therefore, we cannot ensure that the gauge 
coupling in this model, 

A = e*(^), (6.1) 

resembles QCD at the relevant energies. The series solution in Chapter |5] is also 
problematic. This solution is only valid in the regime of z, Zh < 1- However, small 
z corresponds to high energy, the region in which the strong-coupling assumptions 
should fail. 

Our research could be expanded in a number of ways. Since the bulk coordi- 
nate z is related to the renormalized energy scale, one needs to consider running 
parameters as functions of z, 

niq — > mq{z) a — )■ cr(z). (6.2) 

We may even incorporate the anomalous dimensions of operators by assuming 
the scalar mass from the field/operator correspondence is also dependent upon z 
[117] . Further study must be done to sort out the field/operator correspondence 
of the dilaton 0. Whether this scalar field is dual to Tr(F^) or may have 
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consequences to other parts of the theory. In addition, we should investigate if 
the Ajj^ contains relevant information for the gauge theory. The most fertile area 
for further research involves the thermodynamics. We presented a series solution 
that produces many general properties of finite-temperature QCD, but numerical 
analysis may reveal more richness. 

In our research, we attempted to evaluate and improve the relatively simple 
correspondence presented by [H] . Further study into these topics, especially those 
of Chapter |5] will most likely require extensive numerical work and a deeper under- 
standing of stringy effects. While AdS/QCD is not yet a precision tool for QCD 
calculations, it provides qualitative and quantitative insights to difficult problems 
involving strongly coupled gauge theories. 



References 



[1] Particle Data Group, C. A. et al., Physics Letters B667, 1 (2008 and 2009 
partial update for the 2010 edition). 

[2] Y.-Q. Sui, Y.-L. Wu, Z.-F. Xie, and Y.-B. Yang, Phys. Rev. D81, 014024 
(2010), arXiv:0909.3887. 

[3] R. C. Myers and S. E. Vazquez, Class. Quant. Grav. 25, 114008 (2008), 
arXiv: 0804. 2423. 

[4] D. Mateos, Class. Quant. Grav. 24, S713 (2007), arXiv:0709.1523. 

[5] H. J. Kwee and R. F. Lebed, Phys. Rev. D77, 115007 (2008), 
arXiv:0712.1811. 

[6] Jefferson Lab F(pi) Collaboration, V. Tadevosyan et al., Phys. Rev. C75, 
055205 (2007), arXiv:nucl-ex/0607007. 

[7] Jefferson Lab F(pi)-2 Collaboration, T. Horn et al, Phys. Rev. Lett. 97, 
192001 (2006), arXiv:nucl-ex/0607005. 

[8] H. Ackermann et a/., Nucl.Phys. B137, 294 (1978). 

[9] P. Brauel et a/., Phys.Lett. B69, 253 (1977). 

[10] WMAP, G. Hinshaw et al, Astrophys. J. Suppl. 180, 225 (2009), 
arXiv:0803.0732. 



115 



116 

[11] W. de Boer, Prog. Part. Nucl. Phys. 33, 201 (1994), arXiv:hep-ph/9402266. 
[12] J. D. Lykken, (2010), arXiv: 1005. 1676. 

[13] M. Peskin and D. Schroeder, An Introduction to Quantum Field Theory 
(Westview Press, 1995). 

[14] J. "Kapusta, B. Miiller, and J. Rafelski, Quark-Gluon Plasma: Theoretical 
Foundations (Elsevier, 2003). 

[15] R. C. Hwa and X.-N. Wang, Quark-Gluon Plasma 3 (World Scientific Pub- 
lishing, 2004). 

[16] STAR, J. Adams et al, Nucl. Phys. A757, 102 (2005), arXiv:nucl- 
ex/0501009. 

[17] PHENIX, K. Adcox et al, Nucl. Phys. A757, 184 (2005), arXiv:nucl- 
ex/0410003. 

[18] BRAHMS, I. Arsene et al, Nucl. Phys. A757, 1 (2005), arXivrnucl- 
ex/0410020. 

[19] B. B. Back et al., Nucl. Phys. A757, 28 (2005), arXiv:nucl-ex/0410022. 

[20] T. Gherghetta, J. I. Kapusta, and T. M. Kelley, Phys. Rev. D79, 076003 

(2009), arXiv:0902.1998. 

[21] T. M. Kelley, S. P. Bartz, and J. I. Kapusta, Phys.Rev. D83, 016002 (2011), 
arXiv: 1009.3009. 

[22] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri, and Y. Oz, Phys. 
Rept. 323, 183 (2000), arXiv:hep-th/9905111. 

[23] I. R. Klebanov, Int. J. Mod. Phys. A21, 1831 (2006), arXiv:hep-ph/0509087. 

[24] J. D. Edelstein, J. P. Shock, and D. Zoakos, AIP Conf. Proc. 1116, 265 
(2009), arXiv:0901.2534. 



117 

[25] S. S. Gubser and A. Karch, Ann. Rev. Nucl. Part. Sci. 59, 145 (2009), 
arXiv:0901.0935. 

[26] S. S. Gubser, (2011), arXiv: 1103.3636. 

[27] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998), arXiv:hep- 
th/9711200. 

[28] S. S. Gubser, I. R. Klebanov, and A. M. Polyakov, Phys. Lett. B428, 105 
(1998), arXiv:hep-th/9802109. 

[29] E. Witten, Adv. Theor. Math. Phys. 2, 253 (1998), arXiv:hep-th/9802150. 

[30] I. R. Klebanov and E. Witten, NucL Phys. B556, 89 (1999), arXivrhep- 
th/9905104. 

[31] B. Zwiebach, A First Course in String Theory, 1st ed. (Cambridge Univer- 
sity Press, 2005), Cambridge, United Kingdom. 

[32] A. Strominger, Nucl. Phys. B451, 96 (1995), arXiv:hep-th/9504090. 

[33] E. Witten, Nucl. Phys. B460, 335 (1996), arXiv:hep-th/9510135. 

[34] G. T. Horowitz and A. Strominger, Nucl. Phys. B360, 197 (1991). 

[35] T. Gherghetta, p. 263 (2006), arXiv:hep-ph/0601213. 

[36] J. Erdmenger, N. Evans, I. Kirsch, and E. Threlfall, Eur. Phys. J. A35, 81 
(2008), arXiv:071 1.4467. 

[37] P. Breitenlohner and D. Z. Freedman, Phys.Lett. B115, 197 (1982). 

[38] G. 't Hooft, Nucl. Phys. B72, 461 (1974). 

[39] E. E. Jenkins, PoS EFT09, 044 (2009), arXiv:0905.1061. 

[40] E. Witten, Nucl.Phys. B160, 57 (1979). 



118 

[41] J. Polchinski and M. J. Strassler, Phys. Rev. Lett. 88, 031601 (2002), 
arXiv:hep-th/0109174. 

[42] J. Polchinski and M. J. Strassler, JHEP 05, 012 (2003), arXiv:hep- 
th/0209211. 

[43] L. Da Rold and A. Pomarol, JHEP 0601, 157 (2006), arXiv:hep- 
ph/0510268. 

[44] L. Da Rold and A. Pomarol, Nucl. Phys. B721, 79 (2005), arXivrhep- 
ph/0501218. 

[45] J. Erhch, E. Katz, D. T. Son, and M. A. Stephanov, Phys. Rev. Lett. 95, 
261602 (2005), arXiv:hep-ph/0501128. 

[46] C. Csaki, M. Reece, and J. Terning, JHEP 0905, 067 (2009), 
arXiv:0811.3001. 

[47] C. Csaki and M. Reece, JHEP 05, 062 (2007), arXiv:hep-ph/0608266. 

[48] A. Karch, E. Katz, D. T. Son, and M. A. Stephanov, Phys. Rev. D74, 
015005 (2006), arXiv:hep-ph/0602229. 

[49] B. Batell, T. Gherghetta, and D. Sword, Phys. Rev. D78, 116011 (2008), 
arXiv:0808.3977. 

[50] P. Colangelo, F. Dc Fazio, F. Giannuzzi, F. Jugeau, and S. Nicotri, Phys. 
Rev. D78, 055009 (2008), arXiv:0807.1054. 

[51] N. Evans and A. Tedder, Phys. Lett. B642, 546 (2006), arXiv:hep- 
ph/0609112. 

[52] N. Evans, A. Tedder, and T. Waterson, JHEP 01, 058 (2007), arXivrhep- 
ph/0603249. 

[53] F. Zuo, Phys.Rev. D82, 086011 (2010), arXiv:0909.4240. 



119 

[54] A. Zee, Quantum Field Theory in a Nutshell (Princeton University Press, 
2003). 

[55] A. Cherman, T. D. Cohen, and E. S. Werbos, Phys. Rev. C79, 045203 

(2009) , arXiv:0804.1096. 

[56] W. de Paula, T. Prederico, H. Forkel, and M. Beyer, Phys. Rev. D79, 
075019 (2009), arXiv:0806.3830. 

[57] F. Jugeau, (2009), arXiv:0907.5048. 

[58] F. Jugeau, Chin. Phys. C34, 1516 (2010), arXiv:0912.2659. 

[59] M. Huang, Q.-S. Yan, and Y. Yang, Eur.Phys.J. C66, 187 (2010), 
arXiv:0710.0988. 

[60] F. De Fazio, PoS CONFINEMENTS, 128 (2008), arXiv:0812.3218. 

[61] H. Grigoryan and A. Radyushkin, Phys.Rev. D76, 095007 (2007), 
arXiv:0706.1543. 

[62] H. R. Grigoryan and A. V. Radyushkin, Phys.Rev. D78, 115008 (2008), 
arXiv:0808.1243. 

[63] B. Galow, E. Megias, J. Nian, and H. J. Pirner, NucL Phys. B834, 330 

(2010) , arXiv:0911.0627. 

[64] J. Babington, J. Erdmenger, N. J. Evans, Z. Guralnik, and I. Kirsch, Phys. 
Rev. D69, 066007 (2004), arXiv:hep-th/0306018. 

[65] S. Hong, S. Yoon, and M. J. Strassler, (2005), arXiv:hep-ph/0501197. 

[66] R. McNees, R. G. Myers, and A. Sinha, JHEP 11, 056 (2008), 
arXiv:0807.5127. 

[67] G. A. Ballon Bayona, H. Boschi-Filho, N. R. F. Braga, and L. A. Pando Za- 
yas, Phys. Rev. D77, 046002 (2008), arXiv:0705.1529. 



120 

[68] J. Soda, Lect.Notes Phys. 828, 235 (2011), arXiv:1001.1011. 

[69] U. Yakhshiev, H.-C. Kim, and Y. Kim, (2010), arXiv: 1001. 5439. 

[70] H.-C. Kim and Y. Kim, JHEP 10, Oil (2008), arXiv:0808.3639. 

[71] L. H. Ryder, Quantum Field Theory, 2nd ed. (Cambridge University Press, 
1996), Cambridge, United Kingdom. 

[72] A. Krikun, Phys. Rev. D77, 126014 (2008), arXiv:0801.4215. 

[73] D. B. Kaplan, J.-W. Lee, D. T. Son, and M. A. Stephanov, Phys.Rev. D80, 
125005 (2009), arXiv:0905.4752. 

[74] J. B. Kogut and M. A. Stephanov, The Phases of Quantum Chromody- 
namics, 1st ed. (Cambridge University Press, 2004), Cambridge, United 
Kingdom. 

[75] J. D. Bekenstein, Phys. Rev. D 7, 2333 (1973). 

[76] J. D. Bekenstein, Phys. Rev. D 9, 3292 (1974). 

[77] S. W. Hawking, Commun. Math. Phys. 43, 199 (1975). 

[78] T. Springer, (2009), arXiv:0908.1587. 

[79] P. K. Kovtun and A. O. Starinets, Phys. Rev. D72, 086009 (2005), 
arXiv:hep-th/0506184. 

[80] C. P. Herzog, Phys. Rev. Lett. 98, 091601 (2007), arXiv:hep-th/0608151. 

[81] K. Kajantie, T. TahkokaUio, and J.-T. Yee, JHEP 01, 019 (2007), arXivrhep- 
ph/0609254. 

[82] D. T. Son and A. O. Starinets, Ann. Rev. Nud. Part. Sci. 57, 95 (2007), 
arXiv:0704.0240. 



121 

[83] N. Evans and E. Threlfall, Phys. Rev. D78, 105020 (2008), arXiv:0805.0956. 

[84] U. Gursoy, E. Kiritsis, L. Mazzanti, and F. Nitti, JHEP 05, 033 (2009), 
arXiv:0812.0792. 

[85] U. Gursoy, E. Kiritsis, L. Mazzanti, and F. Nitti, Phys. Rev. Lett. 101, 
181601 (2008), arXiv:0804.0899. 

[86] U. Gursoy, E. Kiritsis, G. Michalogiorgakis, and F. Nitti, JHEP 0912, 056 
(2009), arXiv:0906.1890. 

[87] M. Panero, Phys. Rev. Lett. 103, 232001 (2009), arXiv:0907.3719. 

[88] J. Alanen and K. Kajantie, Phys. Rev. D81, 046003 (2010), 
arXiv:0912.4128. 

[89] J. Alanen, K. Kajantie, and K. Tuominen, Phys. Rev. D82, 055024 (2010), 
arXiv: 1003.5499. 

[90] D. Li, S. He, M. Huang, and Q.-S. Yan, (2011), arXiv: 1103.5389. 

[91] S. S. Gubser and A. Nellore, Phys. Rev. D78, 086007 (2008), 
arXiv:0804.0434. 

[92] S. W. Hawking and D. N. Page, Commun. Math. Phys. 87, 577 (1983). 

[93] E. Witten, Adv. Theor. Math. Phys. 2, 505 (1998), arXiv:hep-th/9803131. 

[94] M. Shifman and A. Vainshtein, Phys. Rev. D77, 034002 (2008), 
arXiv:0710.0863. 

[95] B. Batell and T. Gherghetta, Phys. Rev. D78, 026002 (2008), 
arXiv:0801.4383. 

[96] A. Karch, E. Katz, D. T. Son, and M. A. Stephanov, (2010), 
arXiv:1012.4813. 



122 

[97] D. Bugg, Phys.Rept. 397, 257 (2004), arXiv:hep-ex/0412045. 
[98] D. V. Bugg, (2006), arXiv:hep-ph/0603018. 

[99] A. H. Fariborz, Phys. Rev. D74, 054030 (2006), arXiv:hep-ph/0607105. 

[100] D. V. Bugg, AIP Conf. Proc. 1030, 3 (2008), arXiv:0804.3450. 

[101] Y.-Q. Sui, Y.-L. Wu, and Y.-B. Yang, (2010), arXiv: 1012.3518. 

[102] P. Colangelo, F. Giannuzzi, and S. Nicotri, Phys. Rev. D80, 094019 (2009), 
arXiv:0909.1534. 

[103] G. Boyd et ai, Nucl. Phys. B469, 419 (1996), arXiv:hep-lat/9602007. 

[104] D. E. Miller, Phys. Rept. 443, 55 (2007), arXiv:hep-ph/0608234. 

[105] A. Cherman, T. D. Cohen, and A. Nellore, Phys. Rev. D80, 066003 (2009), 
arXiv:0905.0903. 

[106] S. Franco, A. M. Garcia-Garcia, and D. Rodriguez-Gomez, Phys. Rev. D81, 
041901 (2010), arXiv:0911.1354. 

[107] J. I. Kapusta and T. Springer, Phys. Rev. D81, 086009 (2010), 
arXiv:1001.4799. 

[108] D. Vercauteren and H. Verschelde, Phys. Rev. D82, 085026 (2010), 
arXiv: 1007.2789. 

[109] F. V. Gubarcv, L. Stodolsky, and V. I. Zakharov, Phys. Rev. Lett. 86, 2220 
(2001), arXiv:hep-ph/00 10057. 

[110] F. V. Gubarev and V. I. Zakharov, Phys. Lett. B501, 28 (2001), arXivrhep- 
ph/0010096. 

[Ill] D. Dudal, J. A. Gracey, N. Vandersickel, D. Vercauteren, and H. Verschelde, 
Phys. Rev. D80, 065017 (2009), arXiv:0907.0380. 



123 

[112] D. Vercauteren, D. Dudal, J. Gracey, N. Vandersickel, and H. Verschelde, 
Acta Phys. Polon. Supp. 3, 829 (2010), arXiv:1011.3649. 

[113] D. Vercauteren and H. Verschelde, Phys. Lett. B697, 70 (2011), 
arXiv: 1101.5017. 

[114] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl. Phys. B147, 

385 (1979). 

[115] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl. Phys. B147, 
519 (1979). 

[116] M. A. Shifman, A. I. Vainshtein, and V. I. Zakharov, Nucl. Phys. B147, 
448 (1979). 

[117] A. Vega, I. Schmidt, T. Gutsche, and V. E. Lyubovitskij, Phys. Rev. D83, 
036001 (2011), arXiv:1010.2815. 

[118] R. L. Arnowitt, S. Deser, and C. W. Misner, Phys.Rev. 122, 997 (1961). 

[119] S. W. Hawking and G. T. Horowitz, Class. Quant. Grav. 13, 1487 (1996), 
arXiv:gr-qc/9501014. 



Appendix A 



Mathematical Methods 



A.l Useful Relations 

Bessel Functions: 

Jn\x ^ oo) ~ \ — COS a; (A.l) 

1^ / N . ( n-K 'K\ , . ^, 

F^x^oo ~ W— sm a;- — -- A.2 

V TTX \ Z 4/ 

Incomplete Gamma Functions: 

Vis.x^ 0) ~ r(s)-^^ (A.3) 

s 

r(s,x^oo) ~ -x'-^e-'' (A.4) 

r(s,0) = V{s) (A.5) 

A. 1.1 Common Differential Equation 

When solving mass eigenvalue equations in the AdS/QCD soft-wall model, we 
often face a differential equation of the form 

r{z) + [z^ + V'(^) = Ei^i^)- (A.6) 
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The eigenfunctions, ipn, and eigenvalues, E, are found in 



2 , / o-n' 

Mz) = e-'^z"^+-^J——L-{z'), (A.7) 

y [m + ny. 

where L™ are associated Laguerre polynomials. 

A. 2 Schrodinger Transform 

This transform is used to eliminate all first order derivatives from a second-order 
differential equation of the form, 

-^"(z) + uj'{z)^'{z) + W{z)^{z) = m2^(z), (A.8) 

where (') denotes differentiation with respect to z |1H]. One can substitute 

^ = e"^^, (A.9) 
^' = e"^ i-^' + io' ifj) , (A.IO) 



into (A.8) so that it becomes 



r.'2 II 

-ip" + —ip - —ip + W{z)i) = m'^i}. (A.12) 
Thus, we have a Schrodinger-like equation, 

-%p" + V{z)iIj = m^%p, (A.13) 

where 

V{z) = "^-"^ + W{z). (A.14) 
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A. 3 Numerical Methods 



A.3.1 Shooting Method 



Since most eigenvalue problems involving equations similar to (A. 13), where ip — )■ 
ipn and m — )■ m„, cannot be solved analytically, and thus, do not give a closed 
expression for m^, a numerical shooting method was employed. This shooting 
method involves iterating through all values of m„ and numerically solving for 
the wavefunction, ipn- For values of m„ not equal to the eigenvalue, ipn diverges; 
the direction of the divergence is the informative piece of information. As we 
iterate through wavefunction solutions, we discover that if successive values of m„ 
produce ipnS that switch their direction of divergence, then the true eigenvalue 
lies between those values. The excitation mode, n, whose ground state begins at 
n = 1, is then determined by the number of antinodes, a, that exist in the function 
where n = a. Figure A.l illustrates the shooting method concept. 



m < Xi , 

m < Xi 




Xi < m < X2 . 

Figure A.l: The two dotted lines show the wavefunction as far away from any 
eigenvalues Aj. The first dotted line has no antinodes and is below the first eigen- 
value. The second dotted line has one antinode but still diverges; therefore, it is 
between the first and second eigenvalues. The dashed line shows a wavefunction 
for m approaching the first eigenvalue, however, it has not changed its direction 
of divergence. The two solid lines represent eigenfunctions corresponding to the 
ground state (one antinode) and the first excited state (two antinodes). 



127 



A.3.2 Matrix Method 



The equations of motion can be reduced to a set of second-order differential equa- 
tions, 



-(f)" + Viiz)(j) + fiz)7T = 0, 

-n" + V2{z)7r + g{z)^ = 0, 



(A.15) 
(A.16) 



where the eigenvalues are contained within the coefficient functions. These equa- 
tions can be expressed as a system of first-order differential equations 



where W is the matrix 



W 



^' + W{z)^ 


= 0, 




( ° ^ 





o\ 


Vi{z) 













1 


V 9{z) 


V2{Z) 


0/ 



(A.17) 



(A. 18) 



and $ is the vector 



(A.19) 



V -< J 

that forms an orthonormal basis of solutions. We can propagate the solution $ 
between two boundary points 

U{z,Zi,Zo,ml)^{zQ), 



(A.20) 



where we solve (A.17) with the appropriate boundary condition at Zq. The eigen- 
vectors and eigenvalues of U are then calculated. We find two large and two small 
eigenvalues corresponding to two nonrenormalizable and two normalizable eigen- 
functions, respectively. Let us assume the eigenvectors U3 and M4 correspond to 
the small eigenvalues, A3 and A4. Then, any solution for $j can be written as 



aUy, + /9'U4, 



(A.21) 
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where we take the boundary condition as 



$i(2o) 



/ 0(-2o) \ 
-0'(^o) 
7r(2o) 



In order for a and /3 to be nontrivial, we must satisfy 



Un 




for Dirichlet or 




(A.22) 



(A.23) 



(A.24) 



for Neumann boundary conditions. We do this by cychng through eigenvalues 



that minimize the determinant of the 2x2 matrix in (A.23) or (A.24). Practically, 
we find the singular points in the graph of the quantity u^v^—vi^u)^ (or uluf — u'^ul) 



versus m^. An abrupt change in its behavior signals an eigenvalue. Of course, the 



elements chosen from the eigenvectors M3 and U4 are dependent upon the choice 
of Neumann or Dirichlet conditions on the boundary zq. 



Appendix B 



Alternative Parameterization for 
the Scalar VEV 



In this section, we explore two different parameterizations of the scalar VEV. 
While the mass spectra are modified, we see that the results are qualitatively 
the same. The parametrizations that we investigate take the form of a Pade 
approximant and a Gaussian function, 

v3 



aiz + a2Z'- 
1 + a-iZ^ 

VG = biz + b2ze'''\ (B.2) 



B.l Pade Approximant 



Using the form (B.l), we find that the desired limiting behavior of the VEV is 
satisfied, 



vp{z — )■ 0) = aiz + (a2 — 0103)2;^ 
02 



vp{z — )■ 00) 



-z. 



as 



(B.3) 

(B.4) 
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According to the correspondence dictionary, the coefficients take the values 



ai 



02 



03 
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2\/3g5fi — 3mgy/K 

Given this particular form for the scalar VEV, the dilaton becomes 

d(j)p{z) 1603(0103 — a2)z^ + K.{ai + a2Z^)^z 

dz 2(1 + a^z"^) (fli + (3a2 — aia^)z'^ + a2a^z^) ' 

with asymptotic behavior, 



^p{z 0) 



Suml ^ 2(r^gj{3mgK + 2V3^g5fi) 



9{3m^K — Ag^mgiJ? 



9 9 
/xV + .... 
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(B.5) 
(B.6) 

(B.7) 



(B.8) 



(B.9) 

(B.IO) 

(B.ll) 
(B.12) 



In the large-^ limit, the leading behavior of (pp is the quadratic term /i^z^ as 
expected. Using the same equations of motion and methods detailed in Chapter 
|4| we find the scalar, pseudoscalar, vector, and axial-vector meson mass spectra 



produced by (B.l), listed in Table B.l 



B.2 Gaussian Function 



Using the form (B.2), we again find the desired limiting behavior of the VEV, 



vciz 0) 



vciz — )• 00) 



h + b2)z + h2hz'^ 



biz. 



(B.13) 
(B.14) 
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n 


ruf, (MeV) 


(MeV) 


rup (MeV) 


rria, (MeV) 


1 


440 


140 


390 


940 


2 


880 


1194 


840 


1290 


3 


1180 


1499 


1150 


1550 


4 


1430 


1746 


1400 


1770 


5 


1650 


1993 


1620 


1970 


6 


1840 


2255 


1820 


2150 


7 


2020 


2528 


2000 


2310 


8 


2190 


2812 


2170 


2460 



Table B.l: Meson mass spectra produced by the Pade Approximant. 
According to the correspondence dictionary, the coefficients become 



hi = (B.15) 



b2 = ^^-W, (B.16) 
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Given (B.2), the dilaton takes the form 



dz 2(62 + &ie''^^' - 262^3^^) 

whose asymptotic behavior becomes 



(B.18) 



4 Oi + O2 

- ^ z^ +[ — -j^ — + . . . , (B.20) 



45'5 V 3 3v3Kmg/ 

4>g{z^ 00) = ^/ + ... (B.21) 
= fx^z'^ + .... (B.22) 

We recover the quadratic dilaton in the large- 2; limit. Using the equations of 
motion and the methods described in Chapter |4], we find the meson mass spectra 



for the scalar, pseudoscalar, vector, and axial mesons, listed in Table B.2 
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n 


rrifo (MeV) 


(MeV) 


rup (MeV) 


ma, (MeV) 


1 


660 


141 


450 


1080 


2 


1130 


1402 


1020 


1500 


3 


1440 


1779 


1380 


1820 


4 


1680 


2037 


1650 


2060 


5 


1890 


2251 


1870 


2250 


6 


2080 


2478 


2060 


2420 


7 


2250 


2721 


2230 


2570 


8 


2410 


2982 


2400 


2720 



Table B.2: Meson mass spectra produced by the Gaussian form (B.2). 



Appendix C 



Holographic Thermodynamics 



C.l Superpotential Method at Finite Tempera- 
ture 

An alternative solution technique can also be used to solve the five-dimensional 
Einstein's equations. It generalizes the superpotential method considered in [95] 
to the case of finite temperature [H]. The method presents a simple way to 
convert a system of second-order differential equations into a system of first-order 



equations. First, we must take the metric 3.99 and define a new bulk coordinate 
u, where 

a{z)dz = du. (C.l) 

We express a new metric form, 

ds^ = e2^(")(e^(")dr2 + dx') + e-^^^^du\ (C.2) 
where f = e^, a = e^. The Einstein equations become 

A+^-^' + ^x' = 0, (C.3) 

UA" + 3AF - + X^) + e-^V = 0, (C.4) 

4i + F + ^ = 0, (C.5) 
F 
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where Q denotes derivatives with respect to u. A superpotential W can be in- 



troduced to replace the second-order equation (C.3) by two first-order equations. 
The system of equations then becomes 



A 



6 



>V(0,X), X 



dx ' 



F + 



(C.6) 
(C.7) 



dx 



+ o>V(0,x)'-7^>V(0,x)i^ 



+e-^V(0,x) = O. (C.8) 



Using this system, we can arrive at similar solutions for the functions 0, /, and 
V{(f), x) as in Chapter [s] In our research, we did not find that the superpotential 
method held an advantage over our method. 



C.2 The Arnowitt-Deser-Misner (ADM) Energy 

The ADM energy is a useful definition of energy for gravitational systems ap- 
proaching an asymptotic, well-defined metric at the boundary |118l 1119] . To ver- 
ify the thermodynamic relation E = + T S, we will compute the ADM energy 
for our black-hole solution with respect to the thermal AdS solution. Considering 
a time slicing of the five-dimensional metric in ADM form, 

ds'^ = -N'^dt^ + ^ij{dx' - N'dt){dx^ - N^dt), (C.9) 

with i,j = xi, X2, Xs, z. In our case, we have 

iV = a(^)V7M, (C.IO) 
1 = aizf (^dx' + (C.ll) 

The expression for the ADM energy is given by (HH 1119] , 

E = j d^BN (^V^^/^K - . (C.12) 
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The integral is performed over a three-dimensional surface at the bulk boundary 
embedded in the 4D constant time shce Sf, where ■jind is the induced three- 
dimensional metric. The superscript refers to the dimensionality of the extrinsic 



curvature. Using three-dimensional equivalents of (5.6) and (5.8), we obtain 



E = ^2^sA^)vW) [^/ma\S) - AS)^^ , (C.13) 
where we have matched the two solutions at the AdS boundary as specified in 



5.80 Taking the limit 5 — )■ gives the final expression 



E = —^{Uzf„g)-^g{zH)] =J' + TS, (C.14) 



where we have used the fact that 



C.3 Phenomenological Thermodynamics 

As shown in Chapter [5} obtaining an exact solution for the soft-wall model with 
a black hole is challenging. Thus, it is sometimes advantageous to follow the logic 
of [501 llU2j . where we are not concerned with dynamically generating the soft- 
wall geometry. In this case, we allow total freedom within the potential to satisfy 
all the necessary equations of motion. The advantage, of course, is removing an 
added complication; however, we sacrifice any thermodynamic relations useful in 
calculating the transition temperature, namely J-" = J SdT. 

To begin the phenomenological analysis of bhAdS, we start with the same 
action in the Einstein frame, 

s. ^ -i^/^^-^(^-^(^^)^->x)^-mx)) 

+ ^ [ d'xV^K, (C.16) 



using the same metric form, 



a{zY ( -f{z) dt^ + dx^ + 



dz^ 



We assume forms for the dilaton and the metric warping factor a{z) 



a{z) 



R 

—e , 
z 

* .2^2 



11 z - g{zh)z 



and then consider two of the equations of motion 

a'{z)^ a"{z) 



12- 



6 



a{z)^ 
nz) + 3fiz) 



a{z) 
a'jz) 
a{z) 



0. 
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(C.17) 

(C.18) 
(C.19) 

(C.20) 
(C.21) 



Using (C.20), we solve for the tachyon field xiz)i 



X{z) ^2;W2V 



+ z-" 



2^4 



24^2 _ 9VM^ 



(C.22) 



where in ( |C.22 ), the integration constant from the x equation has been chosen so 
that x(^) Xo{z) as 2 — )■ and (C.22) is the expansion about z = 0. Imposing 

/(^h) = 0, 



the black-hole conditions on (C.21), 

/(O) = 1 

we find that 

fiz) = 1 - il . \ { ^ ^ 1 - 



aiw 



(C.23) 



l-f,{zh)z^ + 0{z^). (C.24) 



The integral expression in (|C.24|) can be evaluated using (|C.18|) to give 

1 



V{z) 



27r 



2^2 



-Erf 



/i ev^e 



Erf 



/i 



3£ 



1 - e 



(C.25) 



137 

where Erf is the error function. The other two independent equations of motion 



(5.36) and (5.37) are satisfied by specifying the needed potential, which could be 
temperature dependent. 

After specifying the field behavior, f{z), and a{z), we need to find an accept- 
able form for Q. Unfortunately, in the phenomenological picture, we have little 
guidance. Requiring T{zh) to go to the conformal limit in the small-^h region is 
the simplest constraint that we could impose. To satisfy the conformal limit, we 
find again that the leading order must take the form 

Q{zH)=g^. (C.26) 



as it did in Chapter 5 Using (C.26), we find that the temperature. 



^^''> - An Viz,) - 47rL3 V{z,) ' ^^■^^> 
can be expanded in the small- z/j limit, 

T(..) = — +(^— -— + (C.28) 

The temperature has interesting behavior as z, — > oo. When g = \/2/3, the 
temperature converges to a constant, 

T{z^, ^ oo) = 4- (C.29) 

2 

When g takes any other value, the temperature decays exponentially, T ~ e~^'>, 
so that the theory has no minimum temperature. Thus, we define a critical value 
for g, 

gcrit = (C.30) 
At gcriti we can express the temperature in its exact form. 
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Figure C.l shows the plot of the temperature with and without a condensate 
function. We see that the behavior is similar until we reach laige-Zh- The case 
when g = displays a hnear growth, indicating a second, unstable black-hole 
solution 




Figure C.l: The case of g = and g = gent are plotted. The large- z/^ behavior 
differs with the g = case exhibiting a second black-hole solution. 



The entropy depends strictly on the form of a{z 

^3 
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We expand this expression in terms of small-Zh or large-T, 



As expected, we obtain the conformal behavior in the high temperature hmit, 
s ~ T^. At (?crit, the entropy becomes quite simple, 



Figure C.2 shows two cases of the s/T^ vs T. Besides a more rapid ascent to the 



conformal limit for g = (^crit, each case shows little distinction. 

We also find the speed of sound in the phenomenological model also goes to 
the conformal limit as the temperature increases, 

. ^2 1 822, ^^y^a 2 2 

-^(^) - 3-9^+18^- 



At ^^crit, we obtain an exact expression of v_ 



2 



= r- (C.37) 

{3 + 2fi^zl){{l + e^'^l{V^fiZhETf{fiZh)-l)) 



Figure C.3 shows the plot of the speed of sound versus temperature for (7 = and 
9 = fi'crit- As with the entropy, we see similar behavior of in each case. 



We use the same general expression for the free energy, J^, (5.82), but using 
boundary matching, the cut-off in the thAdS solution becomes 

where 6 is the cut-off point in the bhAdS solution. As in Chapter [5| we find that 
it reduces to the general expression, 

J^= U{zh) + 2VQg{zh). (C.39) 
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Figure C.2: A scaled entropy is plotted for two different values of g. Whether 
we consider the condensate function, there appears to be little difference in the 
behavior of the entropy. 



Using (yfcrit, we find that a transition occurs at Zh ~ 1.07 GeV ^, which corresponds 
to a transition temperature of Tc ~ 308 MeV. We plot two cases of the free energy 



in Figure C.4 The phenomenological model also indicates a condensation function 
is needed to obtain a transition. 
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Figure C.3: The speed of sound is plotted for two different values of g. Each case 
approaches the conformal limit of 1/3. When g = gfcrit, there appears to be a 
steeper ascent to this limit, however. 
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z (1/GeV) 

h 



Figure C.4: The free energy is plotted for two different values of g. We find that 
only when the condensate function is nonzero (f/ 7^ 0), do we recover a transition. 
For g = (yfcrit, we find the transition occurs around Tc ~ 300 MeV. This value is 
high, but consistent with lattice calculations of Nf = QCD. 



